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Preface 


THIS  IS  THE  FIRST  of  a  SGiies  of  publications  prepared  as  notes 
for  a  course  in  Space  Technology,  given  by  the  California 
Institute  of  Technology  (Caltech)  in  cooperation  with  the 
Jet  Propulsion  Laboratory  (JPL)  from  June  19  to  July  31, 
1964.  The  program  was  sponsored  by  the  National  Aero¬ 
nautics  and  Space  Administration  under  Grant  No.  NsG-598 
and  was  taught  by  engineers  from  industry,  from  Caltech, 
and  from  JPL.  It  is  planned  that  the  complete  set  will  con¬ 
sist  of 

Volume  I  Spacecraft  Systems 

by  L.  H.  Abraham,^  douglas  air¬ 
craft  CO.,  INC. 

Volume  II  Spacecraft  Mechanical  Engineering 

by  James  L.  Adams,  jpl 
Volume  III  Spacecraft  Propulsion 

by  F.  E.  Marble,  caltech 
Volume  IV  Spacecraft  Guidance  and  Control 

by  J.  R.  ScHULL,  JPL 
Volume  V  Telecommunications 

by  J.  J.  Stiffler,  jpl 


1  Dr.  Abraham,  the  author  of  this  volume,  died  June  16,  1964.  His  part  of  the 
Caltech  course  was  taught  by  another  Douglas  Aircraft  Co.  engineer,  Mr. 
Duncan  Pitman,  who  also  checked  the  author’s  proofs  for  this  publication. 
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CHAPTER  1 


Historical  Preface 


THE  ASPIRATION  TO  CONQUER  SPACE  is  probably  as  old  as  man  him¬ 
self.  It  can  be  identified  with  his  longing  for  survival  or  for  the 
possession  of  supernatural  powers.  Throughout  the  ages  poets,  phi¬ 
losophers,  and  magicians  have  been  entranced  by  the  notion  of 
attaining  the  supernatural  power  of  levitation.  Writers  of  the  17th 
and  18th  centuries  conceived  space  flight  as  a  form  of  escapism  or 
lesson  in  morality.  Cyrano  de  Bergerac  and  others  in  the  17th  and 
18th  centuries  used  the  concept  of  voyages  into  space  as  philosophical 
and  literary  pretext. 

As  the  early  astronomers  began  to  piece  together  the  true  nature 
of  the  cosmos,  writers,  scientists,  and  philosophers  began  to  envision 
other  planets  in  terms  of  their  own,  inhabited  by  intelligent  beings. 
The  famed  scientist  Kepler  is  credited  with  one  of  the  first  imaginative 
essays  on  space.  His  Somnium  relates  a  voyage  to  the  Moon  or 
“Levania’'  by  his  hero,  Duracotus. 

Similar  tales  were  produced  during  the  Renaissance  and  the  17th 
century.  Milton's  Paradise  Lost  and  Gulliver's  third  voyage  to 
Laputa  both  have  space  travel  themes.  Cyrano  de  Bergerac's  hero 
in  The  States  and  Empires  of  the  Sun  travels  to  the  stars  without  wings, 
by  supernatural  powers  in  an  ingenious  machine  which  might  be  in¬ 
terpreted  as  using  jet  propulsion.  De  Bergerac  also  toyed  with  the 
idea  of  launchings  using  gunpowder  and  fireworks. 

Huygens,  in  his  Cosmotheoros,  in  addition  to  giving  us  a  detailed 
description  of  the  planet  Saturn,  also  discourses  on  its  inhabitants. 

Throughout  the  17th  and  18th  centuries,  tales  of  space  travel 
continued  to  be  expressed  as  imaginary  adventures  either  contrast¬ 
ing  with  or  consistent  with  the  contemporary  world.  But  even  during 
this  period  some  scientific  minds  began  to  contemplate  seriously  the 
possibility  of  space  travel.  Wilkes,  for  example,  produced  a  book 
entitled  The  Discovery  of  a  New  World  with  “a  discourse  concerning 
the  possibility  of  a  passage  thither."  This  book  was  first  published 
in  1640  and  is  considered  a  remarkable  product  of  its  time.  In  one 
passage  Wilkes  theorizes  that  the  Moon's  surface  is  covered  with  a 
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pumice-like  material.  The  reasoning  leading  to  this  conclusion  may 
be  regarded  as  somewhat  devious  by  modern  standards;  perhaps,  it 
was  just  a  lucky  guess,  but  nonetheless  it  was  an  indication  of  the 
author^s  insight.  In  another  section  he  conjectures  on  the  composi¬ 
tion  and  distribution  of  meteoroids  in  the  universe,  not  an  uncommon 
pastime  today.  Wilkes  was  even  curious  enough  to  crush  and  examine 
a  meteorite;  this  was  a  remarkable  step  for  his  day. 

The  following  excerpt  from  this  book  speaks  for  itself: 

Too  bold  was  he,  who  in  a  ship  so  frailc,  first  ventured  on  the  treacherous  waves 
to  sailc.  And  yet  now,  how  easie  a  thing  is  this  even  to  a  timorous  and  cowardly 
nature?  And  questionlesse,  the  invention  of  some  other  means  for  our  con- 
vcince  to  the  moon,  cannot  seem  more  incredible  to  us  than  this  did  at  first  to 
them  and  therefore  we  have  no  just  reason  to  bee  discouraged  in  our  hopes  of  the 
like  succssc. 

Yea,  but  (you  will  say)  there  can  be  no  sailing  thither,  unless  they  were  true 
which  the  poets  does  but  faine  that  she  made  her  bed  in  the  sea.  Wee  have  not 
now  any  Drake  or  Columbus  to  undertake  the  voyage,  or  any  Daedalus  to  invent 
a  conveince  for  the  ayre. 

I  answer,  though  wee  have  not,  yet  why  may  not  for  ceeding  times  rayse  up 
some  spirits  as  eminent  for  new  attempts  and  strange  inventions,  as  any  that 
were  before  them.  Tis  the  opinion  of  Kepler,  that  as  soon  as  the  art  of  flying  is 
found  out,  some  of  their  nation  will  make  one  of  the  first  Colonies,  that  small 
transplant  into  the  other  world  (moon).  I  suppose  his  appropriating  this  pre- 
heminence  to  his  own  countrymen,  may  arise  from  an  over  partiall  affection  to 
them.  But  yet  thus  far  I  agree  with  him  that  whenever  that  art  is  invented  or 
any  other  whereby  a  man  be  conveyed  some  twenty  miles  high,  or  thereabouts, 
then  tis  not  altogether  impropable  that  some  or  other  may  be  sucessefull  in  this 
attempt. 

The  fertile  mind  of  Jules  Verne  produced  the  most  important 
scientific  novels  of  the  19th  century.  His  machines  were  rather 
naive  but  he  intuitively  grasped  the  fundamental  principle  that 
flight  to  the  Moon  would  not  be  made  as  lifting  flight.  Unfortunately, 
he  considered  ballistic  projectiles  the  proper  means  of  a  lunar  mission 
rather  than  rocketry. 

But  through  the  ages,  while  these  men  were  philosophizing  and 
dreaming  of  space  travel,  engineers  were  planning  machines,  and 
scientists  were  developing  principles  and  amassing  knowledge  by 
which  actual  space  travel  could  be  achieved. 

The  origin  of  the  rocket  propulsion  principle,  by  which  all  present 
day  spacecraft  arc  launched,  may  be  traced  to  Hero  of  Alexander’s 
aeolipile.  The  use  of  rockets  or  “fire  arrows”  in  the  11th  century  was 
refen'ed  to  by  Nu  Clung  Tsung  Tao,  and  in  the  early  part  of  the  13th 
century,  the  Chinese  put  rockets  to  serious  use  in  repulsing  the 
Mongols.  In  the  16th  century  one  Wan  Hoo  was  killed  in  a  car 
powered  by  47  solid  rockets.  In  the  late  15th  century,  Leonardo  da 
Vinci  sketched  several  devices  employing  rocket  propulsion,  and 
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around  1600,  Giovanni  de  Fontana  designed  and  built  a  rocket- 
powered  marine  torpedo.  In  the  18th  century  the  Indians  used 
rockets,  probably  of  Chinese  origin,  against  the  British.  These 
rockets,  however,  were  so  primitive  and  unreliable  that  the  develop¬ 
ment  of  rocketry  was  largely  ignored,  while  military  engineers  turned 
to  ballistic-type  projectiles  fired  from  cannons. 

Technological  advances  of  the  early  20th  century  opened  up  new 
possibilities  for  rocket  development.  Among  the  modern  pioneers 
to  realize  the  potential  of  rockets  were  the  Russian,  Konstantin 
Tsiolkovskiy,  the  German,  Hermann  Oberth,  and  our  own  Robert 
Goddard. 

Tsiolkovskiy,  as  early  as  1903,  proposed  a  liquid  rocket  engine  and 
during  his  lifetime  was  instrumental  in  arousing  interest  in  rocketry 
in  his  own  country.  Oberth,  despite  considerable  criticism  and  skep¬ 
ticism  from  fellow  scientists,  established  the  foundation  of  flight 
mechanics  which  predicted  long  range  rocketry  and  space  travel. 
Goddard  launched  the  first  liquid  propellant  rocket  on  March  16, 
1926,  which  date  might  well  be  regarded  as  the  date  of  conception 
of  the  ^ 'space  age.'’ 

The  first  serious  use  of  rockets  as  a  military  weapon  occurred  in 
World  War  II.  The  German  V-2,  placed  in  action  in  1944,  was  too 
late  to  play  a  decisive  part  in  the  outcome  of  that  conflict,  but  it 
awakened  the  world  to  the  possibility  of  the  rocket  as  a  weapon. 
The  pursuit  of  superiority  in  these  weapons  was  to  provide  the  means 
for  launching  space  probes,  then  satellites  and,  eventually,  man 
into  space.  Technology  was  now  adequate  and  mental  climate 
favorable  for  space  travel. 

On  October  4,  1957  the  Soviet  Union  put  into  orbit  around  the 
Earth  the  first  man-made  object.  It  weighed  only  184  pounds,  but 
gave  birth  to  the  "space  age."  The  United  States  followed  within 
a  few  months  by  putting  Explorer  I  into  orbit.  These  simple  instru¬ 
ment  packages  orbiting  in  space  were  soon  to  be  followed  by  more 
complex  experiments  and  eventually  by  men. 

Faced  with  the  reality  of  space  flight,  the  United  States  has  pressed 
military  weapons  such  as  Jupiter,  Thor,  Atlas,  and  Titan  into  service 
as  space  boosters.  The  Saturn  system  presently  undergoing  qualifi¬ 
cation  testing  is  the  first  design  expressly  constructed  as  a  space 
booster.  It  is  hoped  that  this  system  will  place  the  first  American 
astronaut  on  the  Moon. 

The  Moon  landing  is  only  the  first  step  into  space,  however.  For 
the  future,  the  sky  is,  literally,  the  limit. 
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CHAPTER  2 


System  Definition 


j  A  SYSTEM  is  defined  as  any  combination  of  parts,  assemblies,  and  sets 
joined  together  to  perform  a  specific  operational  function  or  functions. 

The  process  of  applying  science  and  technology  to  the  study  and 
planning  of  a  system  so  that  the  relationships  of  various  parts  of  the 
system  and  the  utilization  of  various  subsystems  are  fully  established, 
before  designs  are  committed,  is  known  as  systems  engineering,  j  A 
subsystem  is  a  major  functional  subassembly  or  grouping  of  items  or 
equipment  essential  to  the  operational  completeness  of  a  system-  .  ,  \  / 

SPACECRAFT  SYSTEM  COMPONENTS 

Essentially,  a  space  system  consists  of  the  following  components: 

(1)  Boosters — These  are  the  rockets  used  to  set  the  vehicle  in 
motion  before  other  stages  take  over.  Due  to  the  versatility 
of  the  national  space  effort  a  variety  of  launch  vehicles  is 
available,  each  differing  in  launch  capability.  Table  2.1  gives 
some  of  the  dimensions  and  performance  data  on  booster  sys¬ 
tems  in  the  national  inventory. 

(2)  Transfer  stage — This  stage  is  used  to  transfer  the  vehicle  from 
one  trajectory  to  another.  In  interplanetary  travel  it  is  this 
device  which  transfers  the  vehicle  from  an  orbit  of  the  departure 
planet  to  the  target  planet. 

(3)  Spacecraft — The  spacecraft  is  the  vehicle  required  to  perform 
the  mission  after  injection  into  the  desired  trajectory.  The 
spacecraft  may  be  a  simple  scientific  instrument  package  or  a 
complex  manned  system. 

The  Echo  satellite  is  an  example  of  a  simple  vehicle.  It  is  primarily 
designed  as  a  passive  device  whose  total  mission  is  to  provide  a  sur¬ 
face  that  is  reflective  to  radiofrequency  energy.  Its  total  structure 
consists  of  a  balloon  made  of  Mylar  film,  sandwiched  between  thin 
layers  of  vapor-deposited  aluminum.  These  layers  not  only  rigidize 
the  Mylar,  but  also  provide  a  reflective  surface. 

Apollo,  on  the  other  hand,  is  an  example  of  a  complex  spacecraft, 
the  most  complex  yet  conceived.  This  vehicle  must  provide  trans- 

5 


SYSTEM  DEFINITION 


7 


portation  and  life  support  for  a  crew  to  and  from  a  lunar  landing,  and 
then  return  the  crew  to  Earth. 


DESIGN  SEQUENCE  OF  A  SPACE  SYSTEM 


The  evolution  of  a  space  vehicle  design  from  the  mission  require¬ 
ments  to  the  launching  pad  is  shown  in  figure  2.1.  Variations  of  this 
flow  chart  will  exist  within  various  organizations,  but  the  design 
process  is  generally  similar. 


TRANSPORTATION 

REQUIREMENTS 

1 

FIGURE  2.1 — Vehicle  design  sequence. 
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The  design  and  dev^elopment  of  even  the  smallest  space  vehicle  is 
one  of  the  most  complex  efforts  ever  undertaken  by  man.  It  requires 
the  coordination  and  integration  of  a  multitude  of  scientific  and  engi¬ 
neering  skUls  along  with  the  most  prosaic  of  the  everyday  mechanical 
arts.  An  aerospace  organization  employs  the  skills  and  knowledge 
of  people  ranging  from  artisan  to  zoologist  and  uses  such  diverse 
materials  as  aluminum  and  zirconium. 

The  blending  of  these  skills  and  arts  to  produce  ultimately  a  device 
as  remarkable  as  a  space  vehicle  is  a  considerable  achievement. 

The  development  of  the  vehicle  begins  with  a  definition  of  the  basic 
mission  objectives  and  a  range  of  auxiliary  mission  envelopes  within 
which  the  vehicle  might  be  expected  to  operate.  A  serious  and 
imaginative  approach  at  this  time  will  produce  a  vehicle  with  versa¬ 
tility  and  growth  potential  for  future  missions  as  well  as  providing 
those  for  immediate  use.  With  the  volatility  of  the  space  program 
this  versatility  cannot  be  overemphasized. 

MISSION  REQUIREMENTS 

Design  effort  begins  almost  concurrently  with  the  definition  of  the 
basic  mission  spectrum  and  the  completion  of  the  operations  analysis 
of  the  system.  The  ideal  design  development  consists  of  a  series  of 
overall  vehicle  system  design  studies.  The  first  studies  naturally  can 
be  only  conceptual  in  nature.  In  this  phase,  advance  state-of-the-art 
in  propulsion,  structural  design,  materials  and  manufacturing  arts  is 
evaluated  for  advantage  and  possible  incorporation  into  the  design. 
Decisions  made  at  this  time  can  be  most  decisive  in  the  lifespan  of 
the  system.  An  overconservative  selection  of  a  propellant  or  the 
improper  selection  of  construction  material  will  seriously  limit  the 
effective  life  of  the  vehicle.  In  the  extreme,  the  design  can  die 
aborning. 

The  main  purpose  of  this  sequence  of  preliminary  design  studies 
is  to  estimate  the  major  characteristics  of  the  vehicle  and  to  define  its 
performance  and  operational  characteristics.  At  this  time,  necessary 
supporting  systems  are  delineated  and  coordinated.  The  study  will 
yield  the  following  pertinent  data  from  which  the  balance  of  the 
design  may  evolve: 

(1)  Vehicle  performance 

(2)  Payload 

(3)  Staging  optimization 

(4)  Trajectory  shaping 

(5)  Staging  thrust- to- weight  {TjW)  ratios 

At  this  milestone  in  the  design  study,  major  development  problems 
evolve.  The  required  supporting  research  programs  must  now  be 
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considered  both  as  to  schedule  and  risk  involved  for  the  entire  vehicle 
system. 

STAGING 

A  sonde,  satellite,  space  probe,  or  manned  spacecraft  requires,  for 
the  accomplishment  of  its  mission,  a  specific  quantity  of  kinetic  energy. 
Generally  this  is  expressed  in  terms  of  velocity.  For  example,  the 
velocity  necessary  to  sustain  a  circular  orbit  around  the  Earth  is 


where 

Earth’s  radius  (2.09X10^  ft) 

gravity  at  the  Earth’s  surface  (32.2  ft/sec^) 

altitude  of  satellite 

Thus,  to  maintain  an  orbit  at  500  000  ft  a  velocity  of  25  700  ft/sec  is 
required. 

Since  there  is  no  analytical  method  which  allows  equations  describ¬ 
ing  a  multistage  rocket  to  be  solved  for  the  optimum  mass  ratio  of  each 
stage  (if  the  effects  of  gravity  are  considered) ,  the  solution  must  be  one 
of  iteration.  With  rocket  propulsion,  the  velocity  which  may  be 
obtained  is  dependent  on  the  performance  (specific  impulse)  of  the 
propellant,  the  quantity  of  propellant  aboard  the  vehicle,  and  the 
weight  of  the  vehicle  and  payload.  The  weight  breakdown  of  a 
typical  launch  vehicle  is  roughly: 


% 

Propellant  (fuel) _  88 

Propulsion _  4 

Payload _  b 

Structures _  2 


Defining  propellant  ratio  as  X',  then 

y  weight  of  fuel 
gross  weight 

Thus,  for  a  vehicle  with  a  weight  distribution  as  cited  above,  at  the 
start  of  burning  X'  will  be  0.88.  As  burning  progresses,  this  ratio 
decreases  as  shown  in  table  2.2. 

For  the  case  of  staging,  assume  that  the  fuel  is  divided  between 
two  containers  in  the  ratio  of  7:3.  After  the  first  container  is 
exhausted,  it  is  discarded  with  the  first-stage  powerplant  and  any 
auxiliary  structure.  Under  these  circumstances,  the  payload  will 
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FIGURE  2,3 — Effccf  of  stagins  on  burnout  velocity  of  a  launch  vehicle. 


remain  constant  in  weight,  but  the  propulsion  system  weight  as  well 
as  the  structural  weight  will  decrease.  That  is,  after  staging,  the 
weight  breakdown  might  be  as  follows: 

% 


Propellant  (fuel) _ 75 

Propulsion _  4 

Payload _  17 

Structure _  4 


After  staging  the  mass  fraction  would  be  0,75  instead  of  0.68  as  shown 
in  figure  2.2.  The  effects  of  staging  on  velocity  are  shown  in  figure  2.3. 

Actually,  X'  is  a  measure  of  the  efficiency  of  the  system.  The  higher 
X',  the  more  energy  there  is  available  to  accelerate  the  payload. 
Thus,  if  the  inert  structure  could  be  discarded  as  soon  as  it  is  no 
longer  required,  the  overall  efficiency  of  the  system  would  then 
be  maximum.  Dividing  the  fuel  into  a  number  of  containers  and 
discarding  each  as  it  is  emptied  would  be  ideal.  In  practice,  usually 
only  two  or  three  such  stages  are  realistic. 
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Laiincli  dyiiainics  considerations  usually  require  a  takeoff  thrust- 
to-wcight  {TjW)  ratio  of  about  1.2  to  1.3:1  (a  value  as  low  as  1.1  is 
ac(‘,eptablc  if  the  wind  velocity  at  launch  can  be  restricted  to  a  low 
value).  'This  iininediatcly  presci-ibes  tlie  engine  requirements  for 
the  first,  stage.  Second  stage  TjW  ratio  is  generally  influenced  by 
first,  st.age  burnout  velocity  and  can  vary  from  1  to  2  at  the  lower 
burnout  velocities.  Third  stage  TjW  ratio  should  be  of  the  order 
of  about  1:1  for  a  nuiltimission  vehicle.  Optimum  last  stage  T/IT 
ratio  shifts  widely  a(‘cording  to  injection  conditions,  but,  except  for 
unusual  ojicrations,  performance  is  not  greatl}^  affected  for  a  range 
of  values  from  0.0:1  to  1.2:1.  It  must  be  kept  in  mind,  however, 
that  these  factors  are  not  hard  and  fast  and  require  considerable 
study  to  establish  the  most  favorable  ratios. 


CHAPTER  3 


Sizing  the  Vehicle 


THE  INCREASING  SIZE  OF  SPACE  VEHICLES  and  the  increasing  complexity 
of  their  missions  puts  continuously  greater  pressure  on  the  designer 
to  determine  the  most  efficient  design  for  each  vehicle  proposed. 
Before  the  detail  design  of  the  space  vehicle  can  begin,  certain  overall 
system  studies,  as  sequenced  in  figure  2,1,  must  be  completed.  The 
sizing  and  configuration  study  which  creates  a  vehicle  that  will  per¬ 
form  the  flight  mission  is  a  complex  procedure.  Figure  3.1  illustrates 
this  study  process. 

In  sizing  a  vehicle,  the  choice  of  length  and  diameter  {LjD)  ratio  is 
influenced  by  several  parameters.  In  order  to  arrive  at  a  minimum 
structural  weight  each  of  these  parameters  must  be  taken  into  account 
and  the  overall  optimum  structural  configuration  evolved.  Thus, 
the  influence  of  dome  shape  LjD  ratio,  engine  weight,  interstage  and 
skirt  weight  and  pressurizing  system  weight,  to  mention  the  most 
important  items,  must  be  evaluated. 

Figures  3.2  and  3.3  show  the  effect  of  dome  shape  on  the  weight  of 
the  engine  and  the  interstage  and  skirt  of  a  given  stage.  With  the 
shallower  dome  it  is  to  be  expected  that  a  shorter  interstage  and 
skirt  would  be  feasible.  On  the  other  hand,  increased  engine  weight 
in  the  smaller  diameters  will  offset  most  of  this  advantage.  Then, 
as  shown  in  figure  3.4,  a  minimum  weight  LjD  may  be  determined. 

These  trade-offs  are  extremely  difficult  to  handle  manually;  con¬ 
sequently,  most  of  these  evaluations  are  accomplished  by  digital  com¬ 
puters.  The  principal  purpose  of  these  figures  is  to  emphasize  the 
complexity  of  the  sizing  operation  and  the  futility  of  a  manual  at¬ 
tempt  on  anything  except  the  simplest  design. 

Even  after  these  optimizations  are  accomplished,  there  must  be 
certain  logical  perturbation.  For  example,  from  a  control  standpoint 
it  is  desirable  to  keep  the  diameter  as  large  as  possible.  If  this  can 
be  accomplished  without  undue  weight  penalty  or  performance  loss, 
it  simplifies  the  guidance  and  control  system  to  the  extent  that  phasing 
networks  and  body  mode  accelerometers  can  be  omitted  in  a  very 
rigid  vehicle. 
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FIGURE  3.1 — Flow  chart  for  configuration  study. 

In  general  the  sizing  operation  will  follow  the  flow  chart  shown  in 
figure  3.1  In  this  process  each  configuration  must,  in  turn, 

(1)  Be  sized  to  accomplish  the  mission 

(2)  Be  flown  through  its  mission  profile 

(3)  Have  aerodynamic  flight  loads  and  mass  distribution  data 
determined  on  the  basis  of  standard  gust  wind  shear  and  ground 
wind  criteria 

(4)  Have  critical  loads  determined 

(5)  Have  preliminary  stress  analysis  performed  including  tem¬ 
perature  effects  on  materials 

(6)  Have  weiglit  of  components  computed 

(7)  Utilizing  tlie  newly  determined  mass  distribution,  have  a 
new  time-temperature  history  determined 

(8)  With  this  new  time-temperature  history,  have  structural 
thicknesses  adjusted  and  recycled  until  the  process  converges. 
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FIGURE  3.2 — Effect  of  dome  heisht  on  skirt,  interstage  and  engine  weight. 


Engineering  design  is  an  iterative  decision-making  process.  It 
involves  identification  of  a  need  and  an  examination  of  possible 
solutions  which  are  tested  against  recognized  criteria  of  performance. 
Many  repetitions  of  each  step  are  made  in  a  search  for  acceptable  or, 
hopefully,  optimal  results.  Each  step  in  each  iteration  involves 
choices  and  decisions.  The  more  information  which  can  be  utilized 
in  each  decision,  the  better  the  decision  is  likely  to  be.  This  informa¬ 
tion  is  of  two  kinds:  that  which  is  stored  in  the  designer’s  mind 
(experience)  and  is  readily  accessible  for  use,  and  that  which  must 
be  acquired. 

Usually  two  or  three  promising  configurations  will  result  from  this 
preliminary  design  exercise.  All  of  these  designs  will  have  merit 
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FIGURE  3.3 — EFfects  of  dome  shape  on  booster  weight.  Pressure  — 20-30  Ib/in.^ 


and  only  by  further  study  and  evaluation  will  the  selection  of  a  single 
design  be  possible.  Figure  3.5  shows  an  example  of  the  results  of  a 
second  stage  configuration  study.  Note  the  variation  in  configura¬ 
tion  possibilities  that  provides  adeq'uate  mission  conformity. 

Once  this  selection  is  made,  development  and  research  problems 
may  be  pinpointed.  Engineering  and  production  schedules  are  now 
established  and  the  general  sequence  of  order  given  in  figure  2.1  may 
be  set  in  motion.  The  detail  design  then  developes,  with  the  end 
product  being,  hopefully,  a  successful  operational  vehicle. 

A  complete  multistage  vehicle  configuration  study  can  be  performed 
with  these  iterative  procedures.  Each  stage  is  analyzed  in  turn 
until  all  stages  have  been  evaluated  and  further  iterations  fail  to 
produce  significant  changes  in  any  stage. 

Load  sources  acting  on  a  spacecraft  structure  may  be  categorized 
into  two  types:  static  and  dynamic.  Static  loads  are  of  particular 
use  in  evolving  the  preliminary  structural  design.  For  this  purpose, 
the  structure  is  considered  to  act  as  a  rigid  body.  However,  once 
a  preliminary  design  has  been  achieved,  it  is  possible  and  necessary 
to  compute  the  dynamic  characteristics  of  the  structure.  With  this 
delineation  of  the  dynamic  characteristics,  it  is  now  possible  to  in¬ 
clude  the  efihcts  of  the  dynamic  loads  on  the  system  behavior.  This 
implies  study  of  tlie  interactions  between  the  control  system,  propul¬ 
sion  system,  aerodynamic  effects,  and  the  dynamic  behavior  of  the 
structure. 
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Material  =  aluminum. 


CHAPTER  4 


Aerodynamic,  Powerplant 
and  other  Loads 


THE  COMPONEXTS  OF  A  BOOSTER  OT  Spacecraft  must  be  designed  to 
survive  during  manufacturing  operations,  transportation  to  the  launch 
site,  launch,  mission  operation,  and  reentry.  In  addition,  there  may 
be  other  considerations  such  as  storage  or  erecting  procedures  that  may 
impose  special  loads  or  environments. 

Before  the  engine  is  ignited,  but  while  the  operational  vehicle 
(boosters  and  spacecraft)  is  resting  in  the  launch  complex  awaiting 
takeoff,  ground  winds  can  induce  severe  loads  both  as  a  steady  drag 
load  or  as  a  dynamic  response  in  a  direction  mainly  normal  to  the  wind 
direction.  This  latter  response,  due  to  vortex  shedding,  is  a  low- 
velocity  phenomenon  which  will  occur  at  a  critical  velocity  that  is 
dependent  on  the  natural  frequency  of  the  structure,  Reynolds  number, 
placement  of  protuberances,  tunnels,  etc.  Even  the  proximity  of 
ground  support  equipment  in  the  immediate  vicinity  of  the  craft  can 
influence  the  critical  velocity  and  the  severity  of  the  oscillation. 
Figure  4.1  shows  an  example  of  the  increase  in  bending  moment  over  a 
rigid  body  analysis  resulting  from  vortex  shedding. 

At  engine  ignition  and  release,  longitudinal  transient  loads  are 
generated  which  can  be  rather  severe  not  only  as  they  affect  the  basic 
structure,  but  also  as  they  affect  smaller  components  of  equipment. 
The  magnitude  of  these  effects  is  illustrated  in  figure  4.2.  The  rigid 
body  analysis  would  consider  only  a  steady-state  thrust  force  while, 
actually,  the  interaction  between  the  thrust  transients  and  the  elastic 
characteristics  of  the  structure  would  produce  an  amplification  factor 
of  almost  2. 

Engine  noise  is  most  severe  in  close  proximity  to  the  ground;  the 
noise  generated  by  the  engine  can  be  of  extremely  high  intensity.  This 
sound  pressure  will  affect  the  structure  as  well  as  flight  or  ground 
service  equipment.  Considerable  attenuation  of  the  sound  pressure 
is  often  provided  by  the  structure,  but  in  many  cases  even  this  attenua¬ 
tion  is  not  sufficient  to  prevent  damage  to  delicate  instruments  and 
equipment. 
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FIGURE  4.2 — 3<j  axial  load  at  station  153  vs  time  of  flight.  Boosted  weight  — 
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As  the  spacecraft  leaves  the  ground,  noise  diminishes,  not  only 
because  of  the  attenuation  of  reflected  noise,  but  also  because  of  the 
reduced  density  of  the  air  and  Doppler  attenuation.  It  must  be 
remembered  however,  that  even  at  supersonic  speeds  at  high  altitudes, 
noise  and  vibration  still  can  be  transmitted  by  the  structure  and 
through  fluid  columns. 

In  addition  to  the  engine  noise,  the  following  are  sources  of  vibra¬ 
tion  encountered  during  the  launch  phase  of  the  flight : 

(1)  Engine  ignition  shock 

‘(2)  Thrust  pulsation  (chugging) 

(3)  Turbine  chugging 

(4)  Boundary  layer  noise 

(5)  Fuel  sloshing 

(6)  Control  forces 

(7)  Nonstationary  aerodynamic  forces 

(8)  Unbalance  of  spinning  components 

(9)  Engine  burnout  and  staging  shocks 

During  flight  through  the  transonic  regime  to  the  maximum 
dynamic  pressure  region,  the  thrust  and  various  other  local  steady- 
state  oscillatory  aerod3mamic  loads  become  important.  In  this  phase 
of  the  flight,  the  spacecraft  may  experience  transient  as  well  as  steady 
state  accelerations  resulting  from  atmospheric  winds  and  wind  shears, 
and  from  static  pressure  peaks  which  appear  around  geometric  dis¬ 
continuities.  These  latter  peaks  are  especially  aggravating  at 
transonic  speeds  where  they  may  be  superimposed  on  buffeting  loads. 


TIME  (SEC) 


FIGURE  4.3 — Temperature  history  of  the  transition  section. 
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FIGURE  4.4 — Maximum  fcmperaturcs  of  A-12  trajectorv. 


In  this  same  flight  regime,  consideration  must  be  given  to  the 
control  of  the  vehicle.  This  aspect,  as  well  as  stability  and  the 
coupling  of  the  control  sensors  with  the  vehicle  flexibility  modes, 
must  be  considered  in  making  any  necessary  corrections.  Fuel 
sloshing  can  be  an  additional  aggravation  during  transonic  flight. 

Flutter  of  aerodynamic  surfaces,  such  as  fins,  will  generally  be 
critical  in  the  transonic  maximum-dynamic-pressure  regime.  Thin, 
large,  unsupported  panels  may  also  pose  a  flutter  problem. 

As  the  spacecraft  rises  in  altitude,  the  exhaust  plume  of  the  engine 
will  widen  as  the  atmospheric  pressure  is  decreased.  This  increase 
in  plume  size  and  radiation  can  cause  severe  heating  problems  at  the 
base  of  the  structure. 

Another  serious,  and  often  neglected,  consideration  in  vehicle  design 
is  trajectory  variation.  This  can  result  from  many  sources  but  it  is 
sufficient  to  state  here  that  a  realistic  tolerance  must  be  placed  on  a 
trajectory  and  the  ramifications  of  these  tolerances  must  be  con¬ 
sidered.  Figure  4.3  reflects  the  variation  in  temperature  due  to  the 
flight  trajectory  being  on  the  low  side.  Without  previous  considera¬ 
tion  this  effect  could  have  caused  a  structural  failure.  Figure  4.4 
shows  these  maximum  effects  on  various  sections  of  a  Thor  booster. 


CHAPTER  5 


Material  Selection 


IN  SELECTING  THE  PKOPER  STRUCTURAL  MATERIALS  foT  a  SpaCG  VeMcle 

booster  the  three  most  important  problem  areas  which  must  be  con¬ 
sidered  are: 

(1)  The  environment:  The  natural  and  operational  environment  to 
which  the  material  will  be  exposed  must  be  understood  and  the 
combined  effects  as  well  as  the  individual  effects  on  the  material 
must  be  studied 

(2)  Material  properties:  In  general,  the  physical  and  mechanical 
properties  of  most  materials  are  readily  available.  If,  on  the 
other  hand,  better  materials  are  required  to  cope  with  the 
environment  or  the  special  loadings  of  a  given  mission,  a  better 
understanding  of  both  the  nature  and  origin  of  these  properties 
will  often  yield  a  more  desirable  solution.  This  understanding 
may  go  so  far  as  an  understanding  of  such  details  as  molecular 
interactions  which  affect  the  engineering  properties  of  the  ma¬ 
terial.  For  example,  molecular  dislocations  may  cause  fractur¬ 
ing,  or  radiation  effects  may  cause  breaking  of  long-chain 
molecules  in  elastomers.  Such  mutation  could  result  in 
embrittlement  or  loss  of  elasticity  of  the  material 

(3)  Application  to  design:  The  application  ot  a  material  to  a  design 
will  offer  basic  considerations  whereby  a  thorough  and  basic 
understanding  of  physical  phenomena  could  remove  certain 
physical  limitations  imposed  on  a  material.  Ablative  cooling  is 
a  prime  example  of  such  ingenuity.  Cathodic  protection 
against  corrosion  is  another  example 

ENVIRONMENTAL  EFFECTS  ON  MATERIALS 

The  structure  of  a  booster  or  space  vehicle  represents  a  most  unique 
and  unusual  integration  of  materials  and  environment.  With  the 
possible  exception  of  the  nuclear  powerplant,  no  other  device  can 
approach  this  claim  to  complexity.  The  unique  aspects  of  the  space 
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environment  which  must  be  considered  in  any  material  selection  are : 

(1)  Vacuum 

(2)  Particles 

(3)  Meteoroids 

(4)  Operating  temperatures 

In  addition,  the  following  criteria  have  a  direct  bearing  on  the  selec¬ 
tion  of  a  material : 

(1)  Economics  of  the  complete  vehicle 

(2)  Strength-weight  ratio 

(3)  Experience  in  design  and  fabrication  of  the  material 

(4)  General  state-of-the-art  for  production  and  fabrication 

(5)  Temperature  behavior  characteristics 

(6)  Physical  characteristics  when  in  contact  with  propellants 

(7)  Available  insulations  and  bonding  techniques 

Vacuum  Effects 

The  operation  of  a  space  vehicle  may  involve  exposure  to  the  hard 
vacuum  of  space  for  extended  periods  of  months  or  years  rather  than 
minutes  and  seconds.  The  first  implication  of  this  extended  vacuum 
exposure  is  the  loss  of  material  that  may  result  from  evaporation  or 
sublimation.  In  general,  this  loss  will  be  insignificant  in  its  effect 
on  the  load-carrying  ability  of  metallic  elements.  This  will  be  true 
for  most  metallic  materials  even  at  moderate  temperatures.  This 
material  erosion,  however,  will  not  be  insignificant  where  thin  coatings 
are  used  for  thermal  control  or  other  purposes.  Sublimation  can  be 
a  very  serious  problem  where  very  thin  coatings  are  used  to  control 
the  emissivity  or  absorptivity  characteristics  of  a  surface.  With 
very  thin  coatings,  the  loss  of  a  thickness  of  only  a  few  angstrom  units 
can  change  its  sensitivity  spectrum.  Figure  5.1  shows  the  effects  of 
prolonged  vacuum  exposure  on  several  coatings  used  for  radiators. 
These  tests  were  performed  in  a  vacuum  at  a  temperature  commen¬ 
surate  with  the  operating  temperature  of  a  radiator. 

Friction  and  wear  of  a  moving  surface  is  both  a  mechanical  and 
chemical  process.  Most  systems  are  heavily  dependent  on  physio- 
chemical  barriers  for  the  prevention  of  welding  between  the  micro¬ 
scopic  rouglmesses  of  the  moving  surfaces.  Tliese  barriers  are 
usually  reinforced  by  tlie  assistance  of  a  film  of  fluid  lubricant  which 
is  hydrodynarnically  established  and  maintained. 

In  the  absence  of  air  pressure,  liquid  lubricants  will  generally 
evaporate.  A  low-vapor-pressure  fluid  can  be  used  with  a  moderate 
life  expectancy,  but  where  high  temperatures  and  a  long  service  life 
are  expected,  even  low-vapor-pressure  fluids  are  not  totally  satis¬ 
factory.  For  this  application  solid  lubricants  must  be  substituted. 
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With  a  solid  lubricant,  the  faying  surfaces  are  separated  by  a  low- 
vapor-pressure  solid  film.  Such  films  are  not  self -replenishing,  and 
their  effective  service  lives  are  limited.  Even  ball  or  roller  bearings 
are  subject  to  some  sliding  action  which  will  wear>away  solid  films. 
Once  the  vacuum  environment  of  space  causes  this  protective  lubricat¬ 
ing  film  to  break  down,  the  contact  surfaces  become  “clean,”  and 
surface  particles  begin  to  weld  so  that  the  bearing  soon  becomes 
rough  or  even  freezes. 

Many  dry-rubbing  components  depend  on  absorbed  lubricants  to 
reduce  friction;  even  a  material  such  as  steel  is  lubricated  by  its 
oxide  surface  coating.  Once  this  coating  is  broken  under  vacuum 
conditions  where  it  cannot  be  replaced,  the  frictional  forces  will 
increase  greatly  and  under  certain  pressures  the  surfaces  will  weld. 
This  action  is,  of  course,  accelerated  by  increases  in  temperature 
and/or  pressure.  Thus,  in  the  space  environment,  spacecraft  doors 
can  weld  shut,  relay  contacts  stick,  or  switches  freeze. 

Particle  Radiation  Effects 

While  only  a  few  failures  can  be  attributed  to  the  vacuum  environ¬ 
ment,  there  have  been  several  failures  of  spacecraft  that  can  be 
attributed  directly  to  particle  radiation.  Particle  radiation  in  space 
may  be  divided  into  three  categories:  solar  flare  particles,  cosmic 
radiation,  and  trapped  radiation.  The  primary  galactic  radiation 
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consists  of  positively  charged  particles  of  very  high  energy.  But 
even  for  man — the  most  radiation -sensitive  material  with  the  pos¬ 
sible  exception  of  photographic  plates —this  radiation  may  be  negli¬ 
gible  except  after  many  years  of  exposure.  The  energy  spectrum  of 
the  other  major  constituents  of  space  radiation  is  such,  however,  that 
radiation-induced  changes  in  the  materials  can  be  caused  by  ioniza¬ 
tion  and  excitation  processes.  The  coupling  of  these  effects  with  the 
vacuum  environment  will  produce  effects  not  necessarily  predicted 
by  the  addition  of  effects  of  tests  performed  by  separate  exposures. 

Organic  compounds,  which  include  the  man  in  the  system,  are  very 
susceptible  to  the  inonization  and  excitation  processes.  The  muta¬ 
tion  is  generally  proportional  to  the  total  integrated  dose  received, 
regardless  of  the  type  or  quantity  of  the  radiation.  Thus,  the  solar 
flare  may  present  problems  for  some  of  the  more  sensitive  materials 
but  because  of  the  sporadic  nature  of  these  flares  the  time-averaged 
flux  will  be  relatively  light.  Light  weight  shielding  will  suflice  for 
protection  of  most  materials  from  this  source.  Protection  of  the 
hypersensitive  occupants,  however,  will  require  some  other  form  of 
effective  shielding.  Improved  effectiveness  of  passive  shielding  in 
the  form  of  an  active  magnetic  shield  or  a  heavy  passive  shielding 
in  which  the  occupants  may  shelter  during  an  outburst  may  be 
developed.  Kefinements  in  prediction  of  solar  flares  would  make 
flight  postponement  or  cancellation  possible. 

Within  the  Van  Allen  belts,  the  dose  rate  experienced  due  to 
trapped  radiation  will  be  high  enough  to  affect  any  organic  com¬ 
pounds  used  in  the  spacecraft  structure.  However,  most  organic 
structural  materials  require  large  doses  of  such  radiation  before 
macrochanges  occur,  so  long-exposure  protection  from  this  radiation 
may  be  obtained  with  very  light  shielding. 

However,  it  may  not  be  easy  to  devise  adequate  shielding  for  some 
items.  For  instance,  solar  cells  operate  by  radiation  to  which  they 
are  particularly  sensitive.  Shielding  for  solar  cells  and  for  other 
light-sensitive  or  optical  systems  must  remain  transparent  even 
under  radiation.  Most  glasses  discolor  under  radiation  and  as  yet 
there  is  no  adequate  solution  for  this  problem. 

The  principal  emanations  from  nuclear  fission  are: 

(1)  Alpha  particles 

(2)  Beta  particles 

(3)  Fast  neutrons 

(4)  Gamma  rays 

The  radiations  of  principal  concern  to  the  spacecraft  structure  will  be 
the  fast  neutrons  and  the  gamma  rays.  The  fast  neutrons  inflict 
damage  by  what  is  essentially  a  collision  process.  Gamma  rays, 
on  the  other  hand,  are  of  an  electromagnetic  nature  and  produce 
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damage  by  an  ionization  process.  Damage  by  Gamma  radiation 
causes  impurities  and  interstitials  within  the  construction  material, 
and  may  also  create  atomic  vacancies  and  thermal  spikes. 

Metals  and  ceramics  are  mainly  affected  by  neutron  radiation  which 
tends  to  destroy  their  orderly  structural  arrangement.  Other  ma¬ 
terials  may  be  mutated  by  the  destruction  of  chemical  bonds,  creation 
of  free  radicals,  out-gassing,  cross-linking,  and  polymerization.  Natu¬ 
rally,  these  changes  in  composition  and  structure  must  be  accompanied 
by  dimensional  changes.  Of  all  the  materials,  metals  are  the  least 
susceptible  to  this  type  of  damage  while  organic  compounds,  including 
man,  are  most  susceptible.  Plastics  and  elastomers,  as  a  class,  are  also 
quite  susceptible  to  radiation  damage.  Figure  5.2  illustrates  the 
effect  of  out-gassing  time  on  fatigue  life. 

Meteoroid  Protection 

Certainly  damage  to  a  spacecraft  which  would  endanger  the  occu¬ 
pants  or  prevent  fulfillment  of  the  mission  must  be  avoided.  The 
hazard  presented  by  meteoroid  strikes  is  a  difficult  problem  for  the 
space  system  designer  because  of  the  many  uncertainties  in  the  size, 
composition,  weight,  flux  density,  and  velocity  of  meteoroid  particles. 
Even  the  penetration  mechanisms  are  not  well  understood,  and  simula¬ 
tion  of  the  high  end  of  the  velocity  scale  is  very  difficult. 

Considerable  effort  has  been  devoted  to  systems  designed  for  meteor¬ 
oid  protection.  One  of  the  most  promising  devices  is  the  meteoroid 
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bumper.  This  protection  consists  simply  of  a  thin  sheet  of  material 
that  is  placed  in  front  of  the  primary  structure  and  at  some  distance 
from  it.  In  general,  material  selection  plays  little  part  in  the  protec¬ 
tion  of  structures  from  meteoroid  particles.  The  primary  parameters 
here  arc  those  of  geometry  and  the  physical  properties  of  the  materials 
have  only  small  effects.  Where  the  effects  of  the  damage  must  be 
assessed,  however,  the  problem  is  definitely  one  of  material  character¬ 
istics  such  as  fracture  toughness. 

Many  of  the  physical  properties  of  titanium,  steel,  stainless  steel, 
and  aluminum  alloy  systems  are  reasonably  constant.  For  instance, 
all  titanium  alloys  possess  almost  equal  elastic  moduli,  density, 
coefficient  of  thermal  expansion,  specific  heat,  and  conductivity.  This 
is  also  generally  true  for  aluminum  and  precipitation-hardenable 
stainless  steel  alloys.  As  a  result,  the  designer  can  easily  choose 
between  these  alloy  classes,  if  physical  properties  are  an  important 
consideration  in  tlie  vehicle  design.  The  problem  then  is  reduced  to 
finding  the  best  alloy  in  a  given  base  metal  system. 

The  physical  properties  of  most  materials  are  well  documented. 
The  mechanical  properties  that  affect  the  strength  of  a  material 
strongly  influence  its  selection  for  a  particular  application.  These 
particular  properties  are  defined  by  suitable  mechanical  tests  at  the 
appropriate  operating  temperature  of  the  structure  and  are  generally 
available  in  the  literature.  Tensile  or  compressive  stress-strain  curves 
usually  define  these  members  for  the  designer,  and  accurate  repro¬ 
ducible  test  techniques  to  determine  these  properties  are  well  estab¬ 
lished.  A  typical  full-range  stress-strain  curve  with  the  significant 
properties  derived  is  shown  in  figure  5.3. 

BIAXIAL  BEHAVIOR 

If  the  material  under  consideration  is  homogeneous  and  isotropic, 
the  designer  can  estimate  accurately  the  stress-strain  curve  under 
varying  biaxial  fields  using  the  distortion  energy  theory  (the  von 
Mises  criterion)  as  shown  in  figure  5.4.  However,  those  alloys  which 
exhibit  a  hexagonal  close-packed  crystallographic  structure  are  gen¬ 
erally  not  isotropic  in  sheet  or  flat-rolled  form.  This  anisotropy 
is  known  as  texturing.  The  majority  of  titanium  alloys  exhibit 
hexagonal  close-packed  crystallographic  structure,  and  as  a  result  do 
not-  conform  to  the  von  Mises  criterion.  Thus,  titanium  materials 
wliich  tend  to  have  strong  crystallographic  orientation  show  signifi¬ 
cant  differences  in  biaxial  behavior  when  compared  to  isotropic 
materials.  Figure  5.4  also  shows  a  schematic  plot  of  the  yield  surface 
in  the  biaxial  tension  quadrant.  Note  the  significant  increase  in 
yield  strength  as  the  degree  of  thickness  anisotropy  increases,  higure 
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5.5  compares  the  uniaxial,  predicted  biaxial  and  actual  biaxial  curves 
for  a  pressure  vessel  of  titanium  alloy  (Ti-6A1-4V). 

FRACTURE  TOUGHNESS 

Toughness  (fracture  toughness)  is  a  particular  mechanical  property 
that  defines  the  load-carrying  capacity  of  a  material  in  the  presence  of 
a  flaw. 

A  flaw  is  a  material  defect  such  as  a  crack  which  runs  only  partially 
through  the  thickness  of  the  material  and  usually  results  in  plane 
strain  (triaxial  stress)  at  the  front  of  the  flaw.  Load-carrying  capacity 


FIGURE  5.4 — Yield  surfaces  in  biaxial  tension  for  textured  materials. 
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FIGURE  5.5 — Stress-sfrain  behavior  for  anisotropic  titanium.  (16  in.  diam  test  cylinder 

of  Ti-6AI-4V.) 

is  the  gross  or  net  fracture  stress  that  a  material  will  carry  in  the 
presence  of  a  flaw.  Net  fracture  stress  is  the  maximum  load  di\dded 
hy  the  total  cross  section  area  including  the  flaw  area.  A  critical 
flaw  size  is  the  largest  flaw  a  material  can  withstand  before  its  load- 
carrying  capacity  falls  below  its  ultimate  tensile  strength. 

Tlie  most  common  type  of  flaw  found  in  pressure  vessels  is  a  crack 
in  or  near  a  weld  zone  that  penetrates  only  partially  through  the 
thickness  of  the  material.  It  has  also  been  found  that  data  from  uni¬ 
axial  test  specdinens  witli  shallow  cracks  induced  by  fatigue  correlate 
with  the  results  of  actual  pressing  vessel  burst  tests.  Because  of 
this  correlation,  partial-thickness-crack  specimens  can  be  used  in 
choosing  pressure  vessel  materials. 

These  test  results  also  have  direct  value  for  the  designer  and  the 
nondestructive  test  specialist.  That  is,  this  test  quantitatively  evalu¬ 
ates  materials  in  a  manner  tliat  has  direct  usefulness.  Figure  5.6 
covers  two  candidate  materials  for  a  pressure  vessel.  Material  A  is 
significantly  stronger  than  material  B  but  has  a  critical  flaw  size  that 
is  approximately  0.015  in.  deep.  On  the  other  hand,  material  B,  while 
somewhat  weaker  than  material  A,  has  a  critical  flaw  size  approxi- 
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FIGURE  5.6 — Material  selection  based  on  fracture  toughness. 

mately  0.035  in.  deep.  If  ultimate  tensile  strength  were  the  only 
consideration,  material  A  would  be  the  material  selected.  However, 
if  the  nondestructive  test  specialist  could  not  detect  a  flaw  smaller 
than  0.025  in.  deep  with  the  inspection  methods  at  his  disposal,  any 
structure  designed  from  material  A  would  be  unreliable  unless  the 
maximum  design  stress  were  lowered.  On  the  other  hand,  material 
B  would  be  completely  reliable,  since  flaws  less  than  critical  size  could 
be  detected.  Under  normal  circumstances  material  B  would  be  se¬ 
lected  for  the  pressure  vessel  since  it  possesses  the  maximum  attainable 
strength  and  reliability.  If  the  designer  found  that  material  B  did 
not  provide  structural  efficiency  high  enough  to  meet  the  performance 
requirements,  the  problem  would  have  to  be  resolved  by  the  develop¬ 
ment  of  a  more  advanced  detection  technique  which,  if  such  a  tech¬ 
nique  were  possible,  would  generally  be  more  expensive. 
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THE  DESIGN  OF  A  ‘‘shirt-sleeve’^  CABIN  foi  a  Space  trip  or  for  the 
inhabitants  of  a  space  laboratory  requires  a  range  of  cabin  pressures 
from  a  minimum  of  approximately  7  Ib/in.^  to  a  normal  maximum 
close  to  atmospheric  pressure  (14.7  Ib/in.^).  With  the  cabin  so  pres¬ 
surized  considerable  emphasis  must  be  placed  on  preventing  inad¬ 
vertent  decompression  of  the  crew  compartments.  Compartmenta- 
tion  of  the  working  area  gives  an  added  safety  factor,  but  in  the  final 
analysis  the  reliability  of  the  individual  compartments  must  be  ex¬ 
amined  closely  and  made  individually  safe. 

BASIC  FAILURE  MODES 

All  of  the  possible  modes  of  failure  are  essentially  the  same  differing 
only  in  degree.  The  three  basic  failure  modes  whereby  compartments 
may  become  decompressed  are;  explosive  decompression,  fast  leak, 
and  slow  leak. 

Explosive  decompression  is  initiated  by  a  crack  or  flaw  in  the  pres¬ 
sure-compartment  structure  which,  after  reaching  a  “critical  crack 
length,”  will  become  unstable  and  propagate  rapidly.  Cracks  at,  or 
exceeding,  this  critical  crack  length  become  “fast-cracks”  since  energy 
is  absolved  in  the  spontaneous  running  of  the  crack.  In  a  manned 
compartment  this  mode  of  failure  would  in  all  probability  kill  the 
occupants  instantly,  and  for  all  practical  purposes  it  would  destroy 
the  shell’s  geometry  rendering  the  vehicle  irreparable. 

Before  a  crack  reaches  the  critical  length  a  rapid  leak  would  probably 
occur  and  forewarn  the  occupants  of  impending  danger.  This  pre¬ 
supposes,  however,  that  the  initial  crack  is  not  of  critical  length.  The 
energy  stored  in  the  pressurized  compartment  would  be  released  at  a 
rapid  rate  in  the  event  of  a  fast  leak,  but  not  so  rapid  as  to  cause 
secondary  damage  as  in  the  case  of  explosive  decompression.  The 
compartment  bleed-down,  however,  may  be  so  rapid  that  the  crew 
may  be  unable  to  take  timely  survival  action  such  as  donning  pressure 
suits  or  retreating  to  and  sealing  off  a  nondamaged  compartment. 
Hence,  this  mode  of  failure  could  be  fatal  to  the  entrapped  crew  mem¬ 
bers,  even  though  the  vehicle  could  be  repaired  by  the  surviving  crew. 
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The  slow  leak  is  probably  tlie  only  tolerable  pressure  failure.  Cor¬ 
rective  action  can  be  taken  as  long  as  the  bleed-down  rate  is  rather 
slow  and  location  and  repair  of  the  damage  is  possible.  But  even  so, 
this  mode  of  failure  must  be  considered  serious  enough  to  abort  a 
mission  and  cause  an  economic  or  mission  penalty  ranging  from  a 
moderate  increase  in  operating  cost  to  the  complete  abandonment  of 
the  vehicle. 

The  critical  crack  length  (/)  wliich  a  panel  can  sustain  and  still 
remain  stable,  is  de])ende]it  on  the  ratio  of  the  residual  stress  in  a  panel 
of  width  (w)  after  a  crack  has  appeared  (o-r),  to  the  ultimate  strength 
of  the  parent  material  {(Tuu)’  Tliat  is 


_  W 


f^Ult 
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where  A%  is  an  empirical  material  factor  which  varies  with  the  material 
itself.  For  aluminum  this  factor  is  13. 

POSSIBLE  SOURCES  AND  SIZES  OF  FLAWS  AND  DAMAGE 

In  order  to  completely  define  a  fracture-safe  design,  a  more  thorough 
study  of  the  sources  and  causes  of  cracks  and  imperfections  is  neces¬ 
sary.  Undetectable  manufacturing  defects,  operational  accidents 
and  adverse  environmental  occurrences  are  all  possible  sources  of 
such  flaws. 

In  manufacdairing,  inclusions  of  foreign  substances  can  occur  in 
rolled  sheet  and  plat.e  inaterial.  These  inclusions  are  sometimes 
im])ossiblc  t.o  det.ect  even  with  sonic,  X-ray,  or  other  nondestructive 
test  techniques.  These  UTidetected  flaws  as  well  as  faulty  welded 
joints  can  start  cracks.  Weld  joints  are  usually  minutely  inspected 
and  checked  for  leaks,  but  even  so,  some  flaws  of  this  type  are  not 
detectable.  In  areas  of  'Tigh-flaw  risk’’  the  local  stress  is  generally 
maintained  at  a  much  lower  level  than  the  rest  of  the  structure  to 
minimize  the  efl’ects  of  such  imperfections. 

Other  flaws  such  as  dents,  nicks,  and  scratches  occur  in  fabrication 
and  handling  and  create  a  most  frustrating  problem.  In  general, 
these  stress  risers  arc  burnislicd  or  polished  out  but  the  possibility  of 
latent  damage  always  exists. 

A  pilot  error  or  an  operational  failure  of  the  propulsion  control,  or 
guidance  system,  in  a  docking-feri'y-resupply  craft  can  cause  a  col¬ 
lision  resulting  in  a  puncture  of  a  pressurized  compartment  of  the 
spacecraft.  Most  docking  criteria  specify  a  relative  axial  velocity 
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of  about  2  ft/sec  and  a  lateral  velocity  of  about  1  ft/sec.  However, 
even  these  relatively  low  velocities  can  cause  considerable  collision 
shock  and  any  accidental  contact  could  cause  damage. 

Crew  accidents  may  be  responsible  for  other  damage.  A  tool  that 
slips  could  cause  a  small  puncture  and  leak  and,  depending  on  the 
working  stress  of  the  wall,  even  a  small  puncture  could  precipitate  an 
explosive  decompression.  Failures  of  rotating  machinery,  such  as 
centrifuges  and  gyros,  could  also  cause  punctures  of  pressurized  com¬ 
partments. 

The  primary  factors  governing  critical  crack  growths  are:  pressure, 
compartment  geometry,  material,  and  operating  stress  level.  Nearly 
all  of  the  data  accumulated  so  far  concerning  the  critical  crack  lengths 
of  structural  materials  have  been  arrived  at  empirically.  Available 
data  show  the  critical  crack  length  to  be  somewhat  dependent  on 
form  (i.e.,  geometry)  even  for  identical  materials  and  operating 
stresses. 

For  flat  plates  of  various  materials  under  uniaxial  loads  enough 
data  have  been  generated  to  determine  the  envelope  of  critical  crack 
length  as  a  function  of  operating  stress  or  residual  strength.  Suffi¬ 
cient  data  do  not  exist,  however,  to  define  the  envelope  critical  crack 
length  for  biaxially  loaded  elements.  The  best  available  data  indi¬ 
cate  that  biaxial  stresses  are  the  more  critical. 

Although  information  concerning  fracture  mechanics  is  accumulat¬ 
ing  rapidly,  the  state-of-the-art  of  this  technology  does  not  meet  the 
needs  of  current  design  problems.  In  view  of  the  lack  of  knowledge 
concerning  the  catastrophic  nature  of  failure  due  to  cracking,  it  is 
necessary  to  approach  this  solution  in  a  conservative  manner  to  en¬ 
sure  the  structural  integrity  of  the  vehicle  and  the  safety  of  the  occu¬ 
pants.  Transport  aircraft  represent  the  only  source  of  manned  flight 
operational  experience,  relative  to  structural  failure.  Table  6.1 
summarizes  cabin  design  philosophy  concerning  operating  stress 
levels  for  several  current  and  proven  aircraft.  It  is  true  that  aircraft 
cabins  are  subject  to  severe  cyclical  pressure  loading,  which  is  quite 
different  from  spacecraft  cabin  loading,  but  it  must  also  be  pointed 
out  that  a  transport  aircraft  cabin  can  suffer  a  large  tear  and  ensuing 
rapid  decompression  without  catastrophic  results.  It  must  also  be 
emphasized  that  the  proven  reliability  of  aircraft  cabin  design  is 
attributable  only  in  part  to  the  operating  stress  level.  The  atten¬ 
tion  given  to  detail  in  design  practices  and  provisions  for  a  fail-safe 
structure  are  very  significant  factors. 

The  best  available  information  from  past  fail-safe  design  practice 
and  experimental  data  recommend  an  equivalent  fracture  factor  of 
4.0  for  a  manned  space  cabin.  This  is  believed  to  be  commensurate 
with  the  complexities  and  unknowns  involved  in  a  long-term  space 
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TABLE  6.1— Summary  of  Some  Ratios  of  Ultimate  Strensth  to  Working  Stress  for  Current  Aircraft  Pressure  Cabins'^ 
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mission.  This  equivalent  fracture  factor  represents  the  allowable 
tensile  stress  of  the  pressurized  compartment  structural  material, 
divided  by  the  maximum  permissible  stress  at  the  limit  operating 
pressure. 

A  proof  pressure  of  1.33  X  the  limit  pressure  is  used  for  proof  testing 
of  all  commercial  aircraft.  This  pressure  will  adequately  detect 
leaks  arising  from  structural  deflection  around  seals  and  manufactur¬ 
ing  defects  associated  with  the  installation  of  these  seals.  In  addition, 
a  criterion  of  an  ultimate  factor  of  2.0  on  the  limit  pressure  for  the 
external  walls  of  the  cabin  is  required.  The  significance  of  these 
proof  and  \iltimate  pressures,  however,  is  overshadowed  by  the 
equivalent  fracture  factor  of  4.0  for  fracture-safe  design. 

It  should  be  noted  that  structural  components,  formed  by  pressure, 
but  not  subject  to  explosive  decompression  considerations,  are  de¬ 
signed  for  ultimate  pressure.  Examples  of  this  type  of  structure  are 
frames  and  sandwich  bulkhead  faces. 
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Propellant  Tank  Desisn 


OPTIMUM  DESIGN 

THE  PROPELLANT  TANKS  of  a  boostei’  oT  Spacecraft  are  generally 
designed  to  serve  two  functions:  (1)  as  pressure  vessels  and  (2)  as 
structural  members  of  the  booster  or  spacecraft,  for  which  the  tanks 
must  resist  thrust,  shear,  and  bending  moment.  As  pressure  vessels, 
the  main  stresses  encountered  by  the  tanks  are  hoop  stress 

pR 

^hoop  ~  ^ 

and  axial  stress 

_pR 

cy axial- ^ 

Thus,  for  any  given  pressure,  the  wall  thickness  will  be  directly 
proportional  to  the  pressure.  This  is  shown  as  the  pressure  thickness 
line  in  figure  7.1.  The  stress  due  only  to  thrust  acting  on  the  body  is 
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The  stress  due  to  bending  is 
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Thus,  the  maximum  compression  stress  on  the  vessel  wall  may  be 
wi’itten 
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The  pressure  at  which  all  the  compressive  stresses  are  balanced  out 
(neutral  pressure)  is  then 
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Thus,  for  any  given  thrust,  bending  moment,  and  pressure,  there  is 
a  thickness  requirement  for  the  vessel  wall.  This  thickness  required 
is  shown  as  the  monocoque  compression  thickness  curve  in  figure  7.1, 
An  examination  of  figure  7.1  reveals  the  intersection  of  the  thickness 
required  for  pressure  and  the  thickness  required  for  compression 
stability  at  a  single  value  of  pressure.  If  this  pressure  corresponded 
to  the  operating  pressure  of  the  system,  this  design  would  be  optimum. 
Further  examination  of  figure  7.1  reveals  that,  for  an  operating  pressure 
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less  than  the  pressure  at  the  intersection  of  the  pressure  thickness  and 
stability  thickness  curves,  the  vessel  would  be  operating  at  a  very  low 
hoop  tensile  stress.  This  design  then  would  be  compression  critical. 
Likewise,  for  any  greater  pressure,  the  hoop  stresses  in  the  vessel 
would  be  high  or  the  vessel  would  be  tension  critical.  In  this  case,  an 
improvement  in  the  tensile  strength-weight  ratio  of  the  material 
would  pro\nde  a  more  effective  design. 

It  can  be  readily  observed  that,  depending  on  the  bending  moment 
and  the  internal  pressure,  different  methods  of  pressure  vessel  con¬ 
struction  will  be  required  to  achieve  maximum  efficiency.  As  de¬ 
scribed  earlier,  figure  5.3  shows  the  areas  in  which  various  construction 
techniques  are  most  efficient  for  a  pressurized  cylinder  under  thrust, 
bending,  and  internal  pressure.  For  high  internal  pressures  (such  as 
those  occurring  in  a  solid  propellant  vessel),  a  filament-wound  struc¬ 
ture  with  its  high  tensile  strength-weight  ratio  would  be  most  efficient, 
provided  that  the  bending  moment  is  not  too  high.  If  the  bending 
moment  is  too  high,  an  isotropic  wall  structure  is  preferable.  For 
vessels  of  lower  pressures  or  interstages  where  the  pressure  is  zero, 
truss-core  sandwich  and  honeycomb  structures  are  clearly  superior 
except  in  the  very  high  bending  moment  regime.  These  generaliza¬ 
tions,  as  shown  in  figure  7.2,  are  supported  by  existing  design  practices. 
Of  course,  there  are  other  considerations  which  may  completely  in¬ 
validate  such  an  optimization.  For  example,  sheet  and  stringers  may 
be  preferred  as  components  of  a  compression  member;  however,  other 
design  considerations  such  as  discontinuities  and  attachment  problems 
would  preclude  their  use  in  a  pressure  vessel. 

Figure  7.3  is  a  study  of  the  use  of  materials  for  large-diameter  (18  ft) 
tanks.  When  constant  pressure  lines  are  plotted,  it  is  seen  that  most 
materials,  with  the  possible  exception  of  beryllium,  are  compression 
critical  except  at  relatively  high  pressures.  (The  operating  pressures 
of  most  liquid  propellant  systems  are  of  the  order  of  20  to  70  Ib/in.^) 

It  behooves  the  designer  then,  to  improve  the  buckling  efficiency 
of  the  vessel  over  that  attainable  with  monocoque  construction 
(figure  7.1).  Since  Ejb  is  almost  a  constant  (except  for  beryllium), 
improvement  in  compressive  efficiency  must  be  accomplished  by 
geometric  design.  One  such  technique  is  the  use  of  a  waffle-type 
construction,  which  has  been  used  in  the  Thor  and  Saturn  IV  boosters 
for  the  reasons  mentioned.  Figure  7.4  illustrates  the  variation  of 
the  effectiveness  of  a  waffle  cylinder  over  that  of  a  monocoque  cylinder 
achieved  by  varying  the  waffle  geometry.  In  the  Saturn  IV  vehicle, 
selection  of  a  waffle  pattern  for  the  tank  section  afforded  a  weight 
saving  of  approximately  30  percent  over  a  monocoque  structure, 
accounting  for  the  compression  and  bending  requirements.  As  a 
pressure  vessel,  approximately  85  percent  of  the  waffle  weight  was 


PRESSURE  (I 


PROPELLANT  TANK  DESIGN 


43 


STIFFENER  SKEW  ANGLES 


FIGURE  7.4 — Typical  variation  of  stiffener  effectiveness  with  skew  angle.  Waffle- 
milled  integrally  stiffened  cylinders. 


effective  in  restraining  pressure.  This  is  also  indicated  in  the  example 
shown  in  figure  7,1. 

There  are  three  types  of  possible  failure  to  be  recognized  in  con¬ 
sidering  the  buckling  strength  of  an  integrally  stiffened  cylinder:  the 
genei>i  instability  failure  of  the  cylinder  and  two  local  buckling  modes. 
One  of  these  modes  is  a  buckling  of  the  skin  between  ribs  of  stiffeners 
as  shown  in  figure  7.5.  The  other  local  failure  involves  the  crippling 
or  buckling  of  the  stiffeners.  Either  of  these  local  failures  by  itself 
does  not  necessarily  mean  a  total  collapse  of  the  structure  (see  figure 
7.5).  However,  when  local  buckling  does  occur,  it  reduces  the 
theoretical  overall  buckling  strength  and,  at  worst,  could  cause  a 
general  instability  failure. 

The  critical  buckling  stress  for  flat  or  curved  plates  is  of  the  form 


where 


CfcR — K 


/tsV 

12(1  ~M^)  W 


skin  thickness 
X5=rib  spacing 
ju=the  Poisson  ratio 


FIGURE  7.5 — Insfabilify  failure  of  a  cylinder  due  to  buckling  of  the  skin  between  ribs 

of  stiffeners. 

and  K  is  a  function  of  the  diniension  of  the  plate  and  the  bending 
rigidity  of  tlie  stiffeners. 

Test  results  on  skewed  waffle  patterns  show  that  this  equation 
reduces  to  the  form 
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Similarly,  the  local  crippling  of  the  outstanding  leg  of  the  rib  is 
expressed  as 

where 

thickness  of  the  outstanding  leg 
=  height  of  the  outstanding  leg 

A  satisfactory  value  of  the  buckling  coefficient  Kr  has  been  found 
to  be 

i^,=0.416 


The  stress  in  the  skin  and  rib  elements  can  be  found  by  considering 
the  portion  of  the  total  applied  load  carried  by  each  element  (fig.  7.6). 
This  figure  also  illustrates  how  much  of  the  stiffening  material  in  a 
skewed  waffle  is  effective  in  restraining  pressure.  This  is  why  85 
percent  efficiency  in  the  use  of  stiffening  material  to  restrain  pressure 
was  achieved  in  the  Thor  and  Saturn  designs.  The  total  axial  load 
can  now  be  written  as 


or 


cos  d 


Px=Ps^Pt 

N\ 

COS  B 

cos  a  ^ 


cos  B 


2trtw^(TT  COS  B 


All  modes  of  failure  have  been  treated  separately  up  to  this  point, 
and  the  equations  above  provide  an  analytical  technique.  Hence, 
given  certain  geometrical  dimensions  and  material  properties,  the  criti¬ 
cal  buckling  stress  and  load  can  be  determined.  While  these  equations 
are  useful  in  stress  analysis,  they  are  very  cumbersome  for  the  struc¬ 
tural  designer.  The  use  of  equations  in  this  form  requires  that  the 
designer  combine  critical  stress  equations  and  geometric  proportions 


FIGURE  7.6 — Stresses  in  skewed  waffle. 
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and  that  he  produce  an  optimum  configuration  by  an  iterative  process. 
In  order  to  reduce  the  geometric  parameters,  tlie  designer  must  con¬ 
sider  that  the  dimensional  terms  are  normalized  and  working  design 
charts  such  as  shown  in  figure  7.7  can  be  created. 

PRACTICAL  CONSIDERATIONS 

In  many  cases  the  optimum  configurations  represented  by  the 
envelope  curve  will  prove  to  be  impractical  from  a  design  oi  manu¬ 
facturing  viewjioint.  I^or  example,  the  optimum  rib  thickness  may 
prove  to  be  too  thin  or  the  optimum  rib  height  too  great.  The  rib 
height  is  particuilarly  critical  when  the  walTle  sheet  is  to  receive  a 
subsequent  forming  operation  such  as  rolling  into  a  cylinder.  Cost 
can  also  be  a  serious  consideration.  For  the  taller  ribs,  a  greater 
quantity  of  the  original  plate  stock  will  be  discarded  as  cliips.  By 
limiting  the  rib  thickness-height  ratio  a  compromise  between  weight 
and  practical  limitations  can  be  achieved.  Figure  7.7  shows  the 
deviation  from  optimum  design  dictated  by  manufacturing  require¬ 
ments. 

The  compromise  made  in  this  particular  design  will  result  in  buckling 
stress  on  the  rib  which  would  make  skin  buckling  the  first  failure  mode. 
This  is  still  a  balanced  design,  since  failure  will  still  occur  simul- 


FIGURE  7,7— Design  chaU  of  45°  waffle-stiffened  cylinders. 
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taneousl}’'  in  two  of  the  three  failure  modes.  The  design  may  be  con¬ 
sidered  nonoptimum,  however,  since  all  failure  modes  do  not  occur 
simultaneously. 

Thus,  what  started  out  as  a  rather  precise  determination  of  optimum 
weight  and  geometry  has  now  become  involved  with  another  parameter 
which  is  difficult  to  describe,  namely,  practical  rib  dimensions  which 
are  correlated  with  manufacturing  techniques.  Even  such  items  as 
corner  radii  can  produce  a  significant  effect.  Also,  tolerances  can 
have  a  significant  weight  influence  in  the  manufacture  of  waffle-type 
structures. 

Waffle  structure  may  be  manufactured  by  three  techniques — 
chemical  milling,  machining,  and  forging — although  forging  is  but 
little  used.  Chemical  milling  is  a  process  whereby  selective  areas 
are  exposed  to  chemical  reagents  which  progressively  remove  the 
exposed  metal.  This  process  provides  rather  poor  tolerance  control, 
however;  in  fact,  the  tolerances  of  the  chemical  removal  process  must 
be  added  to  the  normal  tolerances  of  the  original  materials.  The 
undercutting  of  an  outstanding  leg  is  also  a  common  fault,  and  the 
radii  produced  as  the  reagent  erodes  the  material  are  rather  generous. 
This  process  increases  the  nominal  weight  of  the  part  because  of  the 
latitude  of  the  tolerances  involved. 

On  the  other  hand,  the  machining  of  the  waffle  structure  from  flat 
plate  stock  pro\ddes  excellent  tolerances  which  are  consistent  with 
normal  machining  operation.  In  addition,  these  tolerances  are  not 
additive  to  the  original  plate  tolerances.  Small  values  of  r/s,  and  rjw 
consistent  with  good  machine  practice,  are  also  achieved. 

Thus,  unless  the  part  must  be  machined  after  forming  to  a  complex 
contour,  as  in  the  case  of  a  tank  dome,  machining  of  the  waffle  pattern 
is  preferable  to  chemical  milling. 

INSTABILITY  FAILURES 

Even  with  the  most  efficient  waffle  or  stiffened  structure,  however, 
there  are  loading  conditions  and  other  considerations  which  dictate 
the  use  of  a  sandwich  structure  (see  figures  7.1  and  7.2).  In  the  selec¬ 
tion  of  a  configuration  for  a  sandwich  structure,  an  attempt  is  made 
again  to  design  for  simultaneous  local  and  general  instability  failure. 

The  local  instability  criteria  of  a  truss-core  sandwich,  as  described 
in  reference  1,  is 


Kiv^Erf 
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wliere 

thickness  of  the  sjindwich  face 
pitch  of  the  truss 

7)  =  plasticity  reduction  factor 

An  evaluation  of  the  buckling  coefficient  K  for  a  single  and  double 
core  configuration  is  given  in  figures  7.8  and  7.9,  respectively.  These 
plots  define  two  inodes  of  instability :  the  buckling  of  the  faces  between 
the  core,  and  the  biuTling  of  tlie  core  itself.  For  optiinuin  design,  the 
selection  of  the  buckling  coefficient  is  based  on  the  simultaneous  oc¬ 
currence  of  these  two  failure  modes  (figure  7.10  shows  the  excellent 
agreement  bet.ween  test  data  and  the  theoretical  method  proposed  by 
Anderson  ^).  Then,  an  optimum  configuration  may  be  achieved  b}^ 
ecpiating  the  general  insttibility  stress  lor  the  buckling  of  the  cylindei 
to  these  local  instability  stresses. 

Idle  general  instability  of  a  truss-core  sandwich  cylinder  for  single 
corrugation  may  be  expressed  as 

1.7UJEtfbf  tan  $  I  T.  ^ 

\6tf  COS 


FIGURE  7.8— Local  bucklins  coefficient  for  sinsle  truss-core  sandwich  plate. 


*  Andci’sou,  M(‘lviii  8.:  Local  I nstahility  of  the  Elements  of  a  Truss  Core  Sandwich 
Plate.  NASA  TU  R-30. 
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FIGURE  7.10— Titanium  alloy  (Ti-6AI-4V)  roll-welded  sandwich  local  instability 

comprcssure  failure. 


wliere 

AF=runTung  load  (Ib/in.) 

Since  the  efFiciency  is  nondiniensional,  tlie  weight  parameter  may 
be  expressed  nondimensiomdly  by  sliowing  it  to  be  a  function  of 
geometry  and  the  “smeared  thickness”  (equivalent  homogeneous  plate 
thickness,  i), 


weight 


Hi) 


In  order  to  express  the  equations  for  local  and  general  instability 
as  functions  of  NjKR  and  tIR  and  to  eliminate  6/  from  the  stability 


PROPELLANT  TANK  DESIGN 


51 


equations,  we  write 


Solving  for  hf 


N 

(j=—. 


t 


hf=Trtf 


12(l-;u^)  \brJ 


(  K^KgEt  Y 


When  the  expression  for  local  buckling  stress  is  equated  to  that  for 
general  instability  and  the  above  relationship  is  substituted  for  hf, 
the  following  expressions  are  derived  for  single  corrugation 


r 


12/3 


er\r) 


sin  9  cos  9 

(^cose+^fj 


y 


j 


and  double  corrugation 


J 

ER  1 


12/3 -N 


^2/3^^1/3  Q  878 


sin  9  cos  9 


For  simplicit}^,  it  is  assumed  that  the  cylinder  buckling  coeflficient 
is  of  the  classical  value  (6^=  0.606).  The  remaining  efficiency  terms 
are  then  maximized  by  substituting  the  proper  values  of  Kg,  Ka,  and 
tjtf  given  by  Anderson.^ 

These  residts  show  that  maximum  efficiency  is  obtained  when 
0=55°,  /c/^/=0.78,  and  Fls=4.23  for  the  single  corrugation;  and 
0=48°,  ^cA;^=0.70,  and  iiLd=4.58  for  the  double  corrugation. 

When  these  values  are  substituted  in  the  above  equations,  the 
following  relationships  are  attained  for  single  and  double  corrugations, 
respectively. 


N 

ER 


=0.432 


5/3 


and 


N 

ER 


=0.467 


6/3 


These  results  indicate  that  the  double  corrugation  is  slightly  more 
efficient  than  the  single  corrugation. 

For  purposes  of  simplification,  the  optimum  design  equations  given 
above  neglect  tlie  effect  of  core  shear  deformation  and  the  variation 
of  the  buckling  coefficient  from  the  classical. 


2  Ibid. 
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FINAL  OPTIMIZATION  EQUATIONS 


For  the  present,  only  the  effects  of  the  core  shear  stiffness  will  be 
discussed. 

Since,  in  practice,  the  effect  of  core  shear  deflection  is  to  reduce  the 
effective  buckling  coefffcient  by  a  sniall  amount,  an  approximate 
method  ()f  solution  will  be  used. 

Timoshenko  and  Gere  ^  express  the  effect  of  shear  deformation 
upon  the  buckling  of  homogeneous  columns  by  the  equation 


1  + 


llPe 

AG 


wliere 


factor 

area 

(7 shear  modulus 
L=length 


This  equation  may  also  be  expressed  as 

AT, 


N.r- 


1+ 


hG 


where  the  shape  factor  is  assumed  to  be  1,  and  hG  is  tlie  transverse 
shear  stiffness  for  homogeneous  plates. 

The  transverse  shear  stiffness  for  a  corrugated  sandwich  structure 
has  two  values,  one  in  the  direction  of  the  corrugated  core,  and  the 
othci*  transverse  to  these  corrugations.  Since  the  shear  stiffness 
transverse  to  tliese  corrugations  is  mucli  smaller  than  in  the  other 
direction,  only  the  smallest  need  be  considered  liere. 

This  transverse  sliear  stiffness  parameter  has  been  evaluated  foi 
both  single  and  double  corrugation,  assuming  hinged  joints,  and  is 
given  by 


ic/?sin^  0  cos  0 


(single  corrugation) 


and 


2icEsh\^  6  cos  d 

(Tv) 


(double  corrugation) 


3  Timoshenko,  Stephen,  and  J.  Gere:  Theory  of  Elastic  Stability.  McGraw-Hill 
Book  Co.,  Inc.,  2nd  cd.,  1901. 
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Substituting  these  equations  into  the  expression  for  shear  stiffness 
gives 

Ne 

Ncr, = - ^  2>,  (single  corrugation) 

'tcE  sin^  d  cos  $ 

and 


2^c£’sin^  e  cos  B 


(double  corrugation) 


For  simplification,  it  is  desirable  to  express  Ne  in  the  denominators  in 
terms  of  the  buckling  coefficient  C.  For  the  buckling  of  a  cylinder 
where  (Tcr={CE*)IR  this  may  also  be  expressed  as 


N- 


CE,H 

R 


where  t*  is  the  equivalent  buckling  thickness  of  the  sandwich. 

When  the  last  three  equations  are  combined  and  the  final  equation 
is  expressed  as  a  function  of  the  buckling  coefficient,  the  expression 


“  1 

1  <*«■(! 

^Rtc  sin^  B  cos  B 

is  obtained  for  single  corrugation,  and 

c  — _ 

1 

— —  4 

^*F(i-m') 

2Rtc  sin^  B  cos  B 

is  obtained  for  double  corrugation.  These  equations  then  reflect 
the  change  in  the  buckling  coefficient  due  to  the  shear  stiffness  of 
the  core. 

This  analysis  is  based  upon  a  semi-empirical  theory  that  is 
dependent  upon  the  experimental  results  of  monocoque  cylinders. 
Since  the  buckling  coefficient  is  a  function  of  the  thickness  to  radius 
ratio  of  monocoque  cylinders,  this  theory  contends  that  this  ratio 
is  also  applicable  to  sandwich  cylinders.  For  sandwich  construction, 
this  thickness  can  be  shown  to  be  a  function  of  the  equivalent 
“buckling  thickness”  t*. 

The  empirical  monocoque  curve  used  here  bases  the  dependency 
of  the  buckling  coefficient  to  the  thickness-to-radius  ratio  on  results 
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FIGURE  7.11— Classical  coefficient  C=0.606. 


given  by  Gerard  and  Becker.^  An  average  imperfection  parameter 
U  equal  to  0.00015  will  be  used.  This  curve  is  plotted  in  figure  7.11. 

For  purposes  of  evaluating  the  variation  in  the  buckling  coefficient, 
it  is  required  that  the  empirical  curve  from  figure  7.11  be  expressed 
analytically.  By  inspection,  a  hyperbolic  function  of  the  form 
y=[xl{a-\-hx)]-\~c  seems  to  describe  the  curve  adequately. 

When  the  constants  of  tlie  above  equation  are  determined  by  a 
method  of  averages,  the  following  equation  is  obtained: 

^-1X10-= 

C=  — - - r*+o-2'^ 

7.52X  10-3+2.76  j 


Simplifying  the  above  equation  gives 


Ge-0.097 


2.24X10^^+! 

0.307X103^+1 

J 


*  Gerard,  George,  and  Becker,  Herbert;  Handbook  of  Structural  Stability.  Part 
III.  NASA  TN-D163,  Sept.  1959. 
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Substituting  this  equation  into  the  previous  results  gives 


and 


0.367X10^-^+! 

<*<(!-/) 

0.097  (2.24X10®^+!^ 

R  tc  sin^  6  cos  0 

1 

0.367X105^+1 

«*?(!-/) 

0.097(2.24x105  ^+1^ 

^2Rtc  sin^  6  cos  B 

The  aim,  of  course,  is  to  substitute  the  above  equations  for  Cct^ 
and  Ccr^  into  the  original  buckling  equations.  It  is  desirable,  for 
purposes  of  optimization,  to  express  these  equations  in  similar  terms, 
that  is,  to  evaluate  t*  in  terms  of  the  other  parameters  expressed  in 
the  original  buckling  equations.  We  write 

*  1.73t/6/tan0  / 

2^+^  Ve^.cos^"^^ 

'^"^cos  0 

and 

1.7.3^^6^tan0  /  2^,  “ 

"  ,  ,  U  \Qt,cos6'^ 

^"■^cos  e 

Solving  the  equated  buckling  and  instability  equations  on  page 
51  for  bf  and  substituting  into  the  equations  for  t*  and  solving  for 
t*/R  gives  for  single  corrugation 

t*  {i.73wy'^  /kX'"  flY'"  1 

R  {12(1-^^)  tan  w  I  /  tc  I  A 

V  V</ cos  0/ V6ty  / 

and  for  double  corrugation 

t*  (1.73ff)"''^  /KX'XtX  1 

R  {12(1-/)  tan9}'«V<?J  \rJ  (^1  |  (  ^  |  1^ 

For  simplification,  the  buckling  coefficient  Ce  in  the  above  equations 
will  be  approximated  by  assuming  the  value  to  be  0.606  (the  classical 
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coefficient).  A  small  error  is  necessarily  introduced  by  this  simplifi¬ 
cation,  but  in  the  region  of  practical  interest  this  error  is  negligible. 
If  we  let 

a =2  cos  e-\rT 


- 1_1 

^  6^^  cos 

M=0.2>0 


solving  for  (hr^,  we  obtain  the  expression 


(01.g,’”cos9/i'Y^ 
(tan  \R) 


0.097 


,an  d)  ^  ^'^a(3  \R ) _ _  ■  1-b'f  / V 

~si\70K,''^cose y  cose  \rJ  \twne/  t, 

(tan  \R)  / 


for  single  corrugation. 


Similarly,  if 


Y=cos  0+- 


+1 

btf cos  6 


M=o.:io 


solving  for  Ccj^^,  we  obtain 


475K/‘^  cos  e 

(tan0)^^VX  \R/ 


l.097|^ 


Uney^yX  Kjl)  + '  ,  1.29  /iV'Y^ V'’  L 

'^Q0K,''^CQse/tV'^  X  sin^  $  cos  e\Rj  \ta.nej  t. 
(tan  ey^''^y\  \R/  J 


for  double  corrugation.  Therefore,  b}^  substituting  Ccts  Ccr^ 
O,  and  simplifying,  we  obtain 


_iV 

ER  \r) 
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and 

The  efficiency  index  is  described  as  the  terms  inside  the  large  bracket. 
Maximizing  the  efficiency  index  results  in  two  solutions  for  the  core 
geometry  for  these  equations. 

These  equations  are  evaluated  and  plotted  in  figure  7.12  for  those 
core  geometries  showing  maximum  efficiency.  Because  of  the  neg¬ 
ligible  differences  in  the  computed  efficiency  between  0=50°  and  55° 
for  the  single  corrugation,  and  between  0=48°  and  50°  for  the  double 
corrugation,  we  ma}^  arbitrarily  select  0=50°  in  both  cases  for  the 
purpose  of  simplification  and  practical  application.  The  maximum 
efficiency  for  the  single  and  double  corrugation  is  then 

e=50°;f'=0.69;  K,=4.20 
and  '' 

0=50°;  ^=0.74;  Ka=4A0 

respectively. 

To  eliminate  the  necessity  of  evaluating  the  buckling  coefficient 
Ccr^  and  Ccr^  each  time  a  solution  is  required,  we  may  express  the 


FIGURE  7.12 — Core  geometry  for  both  single  and  double  corrugation  showing  maximum 

efficiency. 
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curvoR  shown  in  fi^nirc  7.12  aiialytically.  vSince  a  parabolic  curve 
of  the  form  y-^-ax!’  accural oly  doRcribes  the  curve,  the  resulting' 
equations  may  be  expressed  as 


N_ 

ER 


0.000 


1.  7a 


for  single  corrugation,  and 


N 

m 


dl.OOl 


for  double  corrugation. 

Tn  sandwic.il  design  the  allowable  buckling  stress  frequently  exceeds 
tlie  elastic,  limit  of  the  material.  In  this  case  an  iterative  procedure 
is  ne(‘essary  to  evaluate  correctly  tlie  tangent  modulus  of  elasticity 
at  tlie  linai  condition  of  stress  using  the  above  equations. 

This  iterative  pro(‘.edure  may  be  circumvented  if  the  latter  equations 
are  expressed  in  terms  of  the  loading  index  and  the  material  property. 
By  substitnting  the  relation  Njcf^t  and  by  associating  the  tangent 
modulus  strain  by  the  symbol  A,  where  A=alEt,  the  expression 

/A^V- '' 


is  obtained  foi’  single  corrugation,  and 


A<T^’^^==0.m 


■^XO.74 


for  double  corrugation  with  the  aid  of  geometrical  relationships 
shown  below. 


Single  Corrugation 

Double  Corrugation 

t  -A 

t  —  ^ 

4.31 

«,=0.69f/ 

lAMr 

V3 

!>■ 

CO 

ci 

^fA 

Using  these  equations  and  geometrical  relationships,  the  develop 

iimilt  (71  tJlC  V.*  .....  . —  - 

The  critical  buckling  stress  (a)  and  the  face  and  core  thickness  {tf,  t^) 
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for  an  aluminum  alloy  (7075~-T6)  and  a  titanium  alloy  (Ti-6A1-4V) 
are  shown  in  the  nomogram.  The  tangent  modulus  of  elasticity  for 
the  aluminum  and  titanium  alloys  are  shown  in  figures  7.15  and  7.16, 
respectively.  To  eft'ect  a  solution  for  other  materials,  it  is  only 
necessary  to  compute  the  parameter  and/or  for  the 

particular  material  in  question,  and  plot  the  results  in  the  existing 
nomogram. 
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FIGURE  7.1  3 — Optimum  design  of  truss-core  sandwich  cylinders  with  single  corrugation. 
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WHERE  "A’'  -  _1L 


FIGURE  7,16 — Tangent  modulus  of  elasticity  for  titanium  alloy.^  Annealed  Ti-6AI-4V 
alloy.  Typical  value — compression  at  R.T.  [Sacks,  G,,  and  R.  Ford  Pray  III,  cds.. 
Air  Weapons  Materials  Application  Handbook — Metals  and  Alloys.  Air  Research 
and  Development  Command,  ARDC-TR-59-66.] 
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Example 

Given:  N=S  000  lb /in. 
i?=37.5  inches 

Material:  clad  7075-T6  aluminum  alloy  at  room  temperature 
Find:  The  allowable  buckling  stress  and  geometry  for  single  truss- 
core  sandwich  for  optimum  design 

Solution 

Single  corrugation 


N  3  000 
B 


From  Nomogram  Figure  7.13 


(7^=48000  Ib/in.^ 
^=5.1X10-3 
^,=0.014  in. 
^/=0.020  in. 
^=0.062  in. 
1.16^^ 


T= 


0.325  in. 


CHAPTER  8 


The  Common  Bulkhead 


IN  ORDER  TO  SHORTEN  the  interstage  structure  between  tanks  it  has 
become  increasingly  important  to  provide  a  common  bulkhead  or 
dome  between  the  two  vessels.  With  hypergolic  fluids,  a  twin  nesting 
pair  of  bulkheads  is  necessary  to  avoid  catastrophe  in  the  event  of 
a  leak.  A  dual  wall  will  provide  a  degree  of  insulation  for  cryogenic 
fluids.  A  double-wall  common  bulkhead  designed  of  sandwich 
material  provides  a  safe,  insulative  structure.  Since  such  domes 
are  generally  formed  by  pressure  which  tends  to  buckle  the  dome  in 
compression  the  sandwich  material  will  also  provide  an  efficient 
structure. 

The  analysis  presented  here  will  be  concerned  prhnarily  with  estab¬ 
lishing  the  general  buckling  stability  criteria  due  to  an  external  pres¬ 
sure  force.  Local  instability  failures  such  as  wrinkling  or  dimpling 
will  not  be  considered. 

The  buckling  stability  of  sandwich  spherical  caps  will  be  related 
to  the  buckling  of  monocoque  spherical  caps  by  equivalency,” 
which  is  a  method  of  determining  the  required  configuration  of  sand¬ 
wich  construction  which  has  the  same  buckling  stress  as  a  monocoque 
spherical  cap. 

The  buckling  of  monocoque  spherical  caps  can  be  shown  to  be  a 
function  of  the  nondhnensional  geometric  parameter  X  and  the  non- 
dimensional  bucklmg  parameter  P  which  are  defined  as 

(8,1) 


where : 

i?= spherical  radius,  in. 
thickness  of  shell,  in. 
modulus  of  elasticity,  Ib/in.^ 
o=base  radius  of  shell,  in. 
/n=the  Poisson  ratio 
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The  basic  equations  defining  the  buckling  stability  of  spherical 
caps  as  given  in  equations  (8.1)  and  (8.2)  may  be  related  to  sandwich 
shells  once  the  equivalent  buckling  relationships  are  established. 

To  establish  tlie  equivalence  of  X  in  equation  (8.1)  to  that  of  a 
sandwich  shell,  it  is  necessary  to  derive  the  expression  for  X.  Timo¬ 
shenko  and  Gere  ^  have  shown  that  the  critical  buckling  strength  of 
a  monocoque  cylinder  under  axial  compression  may  be  expressed  as 


=D 


PmV  EU  1 
JiU 


(8.3) 


where 


/>=the  bending  rigidity=MV[l2  (1-m^)] 

R  =  t\\e  radius  of  cylinder 
number  of  half-waves 
Z==  length  of  cylinder 

The  general  buckling  equation  of  plates 
expressed  as 

Ktt^D 


and  shells  may  also  be 


(8.4) 


where  K  is  tlie  nondimensional  buckling  coefficient  dependent  upon 
the  shell  and  plate  geometry. 

Equating  equations  (8.3)  and  (8.4)  and  solving  for  K  gives 


K=mH 


Em 


RWm^TT^ 


(8.5) 


Tlie  buckling  of  cylinders  under  axial  compressive  forces  and  the 
buckling  of  s])herical  caps  under  external  pressure  have  been  shown 
by  previous  investigators  as  having  tlie  same  expressions  when  the  base 
radius  (a)  is  substituted  for  the  length  (L)  in  the  stability  equation. 
Therefore,  by  substituting  a  for  L  in  equation  (8.5)  we  obtain 


Ea^h 


Substituting  D—EII{l—ii^)  into  equation  (8.6)  yields 


K^m‘ 


(8.6) 


(8.7) 


i  Timoshenko,  Stephen,  and  Gere,  J.:  Theory  of  Elastic  Stability.  McGraw- 
Hill  Book  Co.,  Inc.,  2nd  ed.,  1961. 
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Since  the  number  of  half-waves  may  be  assumed  to  be  the  same  for 
both  monocoque  and  sandwicli  c^^linders,  it  is  apparent  that  only  the 
term 


need  be  evaluated  to  establish  an  equivalency, 
equal  to  gives 


Setting  this  term 
(8.8) 


Since  hll—llp^  for  a  condition  of  unit  width,  substituting  this  into 
equation  (8.8)  gives 

w  (1— 


or 


VSp 


(8.9) 

(8.10) 


Since  p=  -pr  for  homogeneous  plates, 
Vl2 


X==[12(l-p^)]^^^ 


(8.11) 


Equation  (8.11),  describing  X,  is  the  same  as  equation  (8.1). 

It  has  been  shown  that  the  relationship  between  the  buckling  of  a 
monocoque  cylinder  and  a  sandwich  cylinder  are  a  direct  function  of 
the  radius  of  gyration  p;  therefore,  to  express  X  in  terms  of  a  sandwich 
spherical  cap,  it  is  only  necessary  to  substitute  the  sandwich  radius  of 
gyration  in  equation  (8.10). 

Assuming  the  face  thicknesses  are  equal  and  isotropic,  have  the 
same  modulus  of  elasticity,  and  the  faces  take  the  membrane  stresses, 
the  radius  of  gyration  is 

_ ic'i'if 

2 


(where  the  moment  of  inerta  of  the  individual  faces  is  neglected). 
Substituting  this  value  into  equation  (8.10)  gives 


X,=^[4(1~m^)F^ 


a 


(8.12) 


Equating  equation  (8.11)  and  equation  (8.12)  and  solving  for  h 
gives 

h=l.73(tc+tf)  (8,13) 
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The  rehxtionsliip  between  monocoqne  and  sandwich  structure  from 
equation  (8.13)  was  set  equivalent  to  t*,  or 

t*=h=^1.7S{tc+tf) 

Substituting  in  equation  (8.11)  gives 


X=[12(1-m^)] 


1/4 


a 


(8.14) 


Equation  (8.12)  and  equation  (8.14),  defining  X,  are  equal  for  the 


conditions  assumed. 

To  account  for  face  materials  and  face  thicknesses  which  are  not 
equal  (neglecting  tlie  bending  stiffness  of  the  individual  faces)  the 
more  general  derivation  of  the  equivalent  buckling  thickness  t*  is 
expressed  as 


i*= 


1.73 


'y/ 


+  £*2^ ^2) 


m-uy 


El  U 


(8.15) 


To  evaluate  the  nondimensional  buckling  parameter  P  in  equation 
(8.2)  to  an  equivalent  sandwich  spherical  cap,  the  similarity  of  the 
classical  buckling  equations  for  cylinders  under  axial  compression 
and  spherical  caps  under  external  pressure  will  be  applied. 

According  to  Timoshenko  and  Gere,^  the  buckling  of  monocoque 
spherical  caps  and  cylinders  may  be  expressed  as 

<7„=CE^  (8.16) 


or,  in  terms  of  the  buckling  pressure,  as 


qR 

2h 


(8.17) 


Equating  equations  (8.16)  and  (8.17)  and  solving  for  qcr  gives 

<1„=2Ce{^  (8.18) 

Substituting  equation  (8.18)  into  equation  (8.2)  for  5  gives 

2(1— 


^  Ihid. 
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The  term  —  ij^)\l¥R’^  in  equation  (8.19)  may  be  expressed 
also  in  terms  of  X  by  a  modification  of  equation  (8.8)  as 

h^R‘^  12 


and  substituting  in  equation  (8.19)  yields 


(8.20) 


Since  X  has  already  been  evaluated  for  sandwich  construction  from 
equation  (8.14)  it  is  only  necessary  to  evaluate  the  nondimensional 
parameter  C  in  equation  (8.20)  for  sandwich  construction. 

The  critical  buckling  stress  for  sandwich  cylinders  has  been  shown 
to  be 

(8.21) 


or,  in  terms  of  the  buckling  pressure,  for  sandwich  faces  of  equal 
thickness  and  modulus  of  elasticity 


(8.22) 


Equating  equation  (8.21)  and  equation  (8.22)  and  solving  for  Ogives 


4:Eit% 


(8.23) 


Substituting  C  from  equation  (8.23)  and  X  from  equation  (8.14)  into 
equation  (8.20)  results  in 


(8.24) 


Equation  (8.24)  is  the  equivalent  expression  for  sandwich  construction 
when  the  face  thicknesses  are  equal  and  the  face  moduli  of  elasticity 
are  the  same. 

Equation  (8.24)  expressed  in  terms  of  the  buckling  stress  results  in 


2(l-ju^)aV 


(8.25) 


To  express  equations  (8.24)  and  (8.25)  for  sandwich  structure  of 
different  face  moduli  and  unequal  face  thickness  where  Ei9^E2  and 
ti^t2,  it  is  only  necessary  to  equate  the  buckling  strains  and  satisfy 
the  condition  of  equilibrium. 
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Equating  €1  =  62,  where 


61 


0^1  _  cr2 

E,’ 


gives 


0-2 


O-iE'2^ 


(8.26) 


From  the  condition  of  equilibrium  A^i+A^2=A^  where 

A/"2  =  0'2^2 


gives 


O’!  -|- (72^2  —  N 


(8.27) 


Substituting  equation  (8.26)  into  equation  (8.27)  and  solving,  result  in: 

E2,N 

^  {EiJ,x-\-E2^t2) 

and 

E,^N 

«2) 

Since  N=(jcrEI2^  substituting  this  value  in  the  above  equations  gives: 


and 


0-2  ([ctE 

2{Ei^ti-\-E2^t2) 


0-1  ^  q,rR 

Ei^  2{Ei^ti^E2^t2) 


Substitution  into  equation  (8.25)  results  in 

p_  (1— 

t*'\Ei^ti-^E2j^2) 


(8.28) 


Or  equation  (8.28)  in  terms  of  ci,  o-o  gives 

g'2  29) 

~  t^^RE,^,E2^ 

The  buckling  equations  (8.28)  and  (8.29)  do  not  consider  the  effect 
of  shear  deformation  of  the  core.  In  practice,  the  shear  deformation 
is  of  significance  because  of  the  lack  of  shear  rigidit}^  of  the  core.  To 
account  for  this  shear  deflection,  use  is  made  of  the  discussion  by 
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Timoshenko  and  Gere  ^  of  the  effect  of  the  shearing  deformation  in 
the  buckling  of  solid  isotropic  columns  by  the  equation 


Fe 


i+— 


where 


t^EI 


Fe-  ^2 

shape  factor 
area 

(7= shear  modulus 

Equation  (8.30)  may  also  be  expressed  as 


N.. 


Ne 


^tO 


(8.30) 


(8.31) 


where  the  shape  factor  is  assumed  to  be  1,  and  the  width  1. 

Since  the  buckling  of  a  spherical  cap  and  a  column  are  related, 
equation  (8.31)  will  be  used  to  express  the  effect  of  the  shearing  defor¬ 
mation  upon  the  buckling  of  a  sandwich  spherical  cap.  From  equation 
(8.28)  where 


the  parameter  qcr  will  be  expressed  in  terms  of  Ne  by  substituting 
qcr=‘^NelR  and  solving  for  N.  Therefore, 


Ne 


PF  R{EiJ.i^E2^t2) 


(8.32) 


Letting  the  denominator  l-^NejtG  m  equation  (8.31)  be  equal  to  77, 
equation  (8.31)  may  be  written  as 


The  shear  stiffness  parameter  for  sandwich  construction  equivalent 
to  the  value  of  {tG)  in  equation  (8.31)  is 


tjtPQGc 

2L 


(8.33) 


^Ihid. 
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where  the  sheiir  stiffness  of  the  face  is  neglected.  Substituting  this 
e(|nivalent  sandwich  vahie  of  shear  stiffness  and  Nc  from  equation 
(8.32)  and  (8.33)  and  solving  for  rj  gives 


Pt*’‘RUE:,tt  +  E2^t2) 


(8.34) 


Since  the  faces  may  be  considered  thin  with  respect  to  the  core  thick¬ 
ness  in  most  cases  of  sandwich  design,  it  may  be  assumed  that 

tc^  t 


Therefore,  equation  (8.34)  reduces  to 


77  =  1+- 


2{i-fi^)aHGc 


(8.35) 


From  equation  (8.14)  where  X=[12(l-iu')]'/'  a/^/Rt*,  this  may  also 
be  expressed  as 

Rt*~  12 


a'‘(l  — )u^)  yR 

Substituting  into  equation  (8.35)  gives 


?7  =  l-f 


yRtGc 


(8.36) 


(8.37) 


From  equation  (8.13)  where  ^*=  1.73(^r-|-^/)  assume  that 
Using  tliis  simplifying  assumption,  equation  (8.37)  reduces  to 


77  =  14 


12X1.7374ii’q^i  +  ^2,^2) 


y2GrR 


(8.38) 


Since  Ncr^^Nelv  nuiy  also  be  exjmessed  in  terms  of  the  monocoque 
buckling  coefTicient  substituting  the  value  of  77  from  equation  (8.'38) 
results  in 


Pe _ 

20.7P{E:U-^E2,t2) 

y2GcR 


(8.39) 


Equation  (8.39)  has  been  plotted  in  figure  8.1  to  show  the  variation  of 
the  buckling  parameter  P  as  a  function  of  the  shear  stiffness  of  the 
core.  The  basic  monocoque  curve  was  based  upon  data  using  a  90 
percent  probability  curve. 
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FIGURE  8.1 — Buckling  of  sandwich  spherical  caps  under  external  pressure. 


Example  Problem 

Spherical  radius,  R=S0  in. 

Base  radius,  a-=100  in. 

Face  material,  2014-T6  aluminum  alloy 
Core  material,  HRP  3/16—4  Ib/ft^ 
^c~2.00  in. 

^1  =  0.020  in. 

/2=0.040  in. 

Gc=l2  000  lbs/in.2 
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Outer  face  U  is  at  —423°  F 

Inner  face  U  is  at  room  temperature 

m=0.30 

Neglecting  thermal  gradient  stress,  find  the  critical  buckling  pressure 
qcr. 

From  equation  (8.15)  compute  t*,  assuming  elastic  modulus  values 
of  =12.3X10®  and  7^2  =10.5X10®,  whore 

‘s  ^ 


/ 

1.73 

{Eyfi-]rE2t\)Eiti 

E,^  t,  +  E2^  U  -  [Eyfi  {2t-u)  fiE2t2^f 


Substituting  the  assumed  clastic  secant  moduli  for  Ei^  and  £'2^  at 
temperature  from  figure  8.2,  and  the  value  of  ti,  t2,  and  t  into  the  above 
equation  results  in  ^*  =  3.39. 

From  equation  (8.14)  compute  the  parameter  X,  where 


X=[12(1-m^)] 


1/4 


a 


^/Rt* 


This  results  in  X=  11.1. 

(Compute  the  parameter  {E,fifiE2fi)l{2RGc)  from  figure  8.1,  where 
Ei^  and  E2^  are  the  elastic  tangent  modulus  values.  This  gives 

12.3X10^X. 02  +  10. 5X1Q‘^X.04_ 

2X80X12X10^ 

From  the  curve  in  figure  8.2  reading  in  the  appropriate  computed 
values  gives  Pcr  =  8G0. 

wSiuce  the  buckling  parameter  P^r  uiay  yield  membrane  stresses 
above  the  ])roportional  limit  of  the  facing  material,  it  is  recommended 
that  tlie  buckling  stresses  be  computed  before  the  critical  buckling 
pressure. 

From  equation  (8.29) 


p  _2(1— +)aV2,  0-1 


or  this  may  be  written  as 
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Various  A  values  have  been  constructed  in  figure  8.2.  Substituting 
the  values  of  Per,  t*,  R,  and  a  into  the  above  equation  results  in 

360X  (3.39)^80 
^“2(1-0.3'‘)100* 

ri.=6.17X10-' 

This  value  of  A  described  in  figure  8.2  shows  that  the  predicted  stress 
is  above  the  proportional  limit  of  the  facing  materials;  therefore,  an 
iterative  approach  must  be  used  to  arrive  at  a  solution.  At  A  =  6. 17 
XIO"’,  the  stresses  in  the  faces  from  figure  7.14  are 

<ri=61  500  psi  at  £1, =10X10® 

and 

0-2=51  000  psi  at  £'2, =8.25X10-® 

The  appropriate  secant  moduli  at  these  stresses  are 
£1  =12.1X10®,  £2  =10.3X10® 

Substituting  the  above  values  of  and  E2s  into  equation  (8.I0)  for 
t*  results  in 

i*=3.39 


which  is  the  same  value  as  originally  calculated. 

Since  the  value  of  t*  is  the  same,  and  from  equation  (8.14)  the  value 
of  X  is  a  function  of  the  variable  t*,  therefore  X=ll.l. 

Recomputing  the  parameter  {Eifi-\-E2f2)li2RGc)  using  — 10X10 

and  £’2^==8-25X10®  from  the  first  calculation  results  in  the  value  of 

0.276.  From  figure  8.1,  Pcr=360. 

Since  P  and  Per  remain  the  same,  the  computed  value  of  A  is  also 
the  same.  Therefore,  for  this  particular  example,  a  solution  has  been 
determined  that  satisfies  the  original  values  of  the  secant  moduli. 

To  establish  the  budding  pressure  it  must  be  stated  that,  unlike  a 
solution  in  the  elastic  region  where  both  faces  presumably  fail  to¬ 
gether,  the  plastic  region  results  in  one  face  with  the  minimum  strain 
failing  first.  Therefore,  to  determine  the  actual  failing  pressure,  it 
is  necessary  to  determine  this  minimum  strain.  This  can  best  be  ac¬ 
complished  by  the  following  method 

At  the  A  value  of  6.17X  10-^  the  corresponding  B  values  or  strains 
result  in  P2<Pi-  Since  the  strains  must  be  equal  in  the  faces,  it 
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follows  that  a  solution  exists  at  the  minimum  strain  value  at  B2' 
Therefore  constructing  the  B2  line  in  figure  8.2  where 


B2 


<J2  51  000 

^“10.3X10® 


4.95X10"^ 


and  at  the  intersection  of  the  outer  face  at  —423°  F,  results  in 


0-1=60  000  lb/in.2 

The  allowable  budding  pressure  is,  therefore, 


or 


2=^[60X  10^X0.02+51X10^X0.04] 


2=81  Ib/in.^ 


CHAPTER  9 


The  Tensile  Paradox 


IN  CERTAIN  ELLIPTICAL  DOMES  or  Segments  of  a  truncated  sphere, 
the  application  of  internal  pressure  will  paradoxically  produce  a 
compressive  hoop  stress  which  may  result  in  a  compressive  buckling 
failure  of  the  shell,  even  though  the  shell  is  primarily  under  tension. 
Figure  9.1  shows  a  compressive  budding  failure  on  a  truncated  seg¬ 
ment  of  a  sphere.  A  uniform  tension  load  was  applied  at  the  periphery 
of  the  truncated  segment.  Because  theory  for  this  latter  failure  is 
not  well  represented  in  the  literature,  a  solution  will  be  offered  here. 

A  theory  of  local  instability  is  developed  along  the  lines  of  the  book, 
Non-Linear  Theory  of  Thin  Elastic  Shells,^  edited  by  Mushtari  and 
Galimov,  from  the  approximate  equations  of  compatibility  and  equi¬ 
librium  (equations  57.5  and  57.6  in  the  book) : 

(mi  sin^  9 


DAAw- 


4'ee' 


Wm- 


w  sin 


(9.2) 


It  will  be  shown  at  a  later  stage  of  this  analysis  that  the  assumption 
of  a  large  number  of  buckles  can  be  predicted  by  the  theory  and 
substantiated  by  the  test  results.  In  fact,  the  actual  buckled  con¬ 
figuration  has  been  used  as  a  guide  to  reduce  the  number  of  inde¬ 
pendent  variables  from  two  to  one.  The  number  of  buckles  will  be 
computed  to  ensure  the  validity  of  the  equations  and  this  number 
will  be  compared  with  the  actual  number  as  determined  by  test.  The 
spherical  shell  segment  as  shown  is  loaded  by  tangential  tensile  loads 
at  points  a  and  h  which  are  radially  restrained.  A  single  buckle 
forms  between  a  and  h  and  is  repeated  a  considerable  number  of  times 
around  the  circumference. 


^  Mushtari,  Kh.  M.,  and  Galimov,  K.  Z.,  editors:  Non-Linear  Theory  of  Thin 
Elastic  Shells,  Academy  of  Sciences,  USSR,  1957. 
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FIGURE  9.1 — Observed  buckle  confisuraHon. 


In  order  to  pass  through  points  a  and  6,  the  radial  deflection  function 
^y=C'sin  m(0— ^o)  cosn^ 


must  be  chosen  such  that  m(7r— 2^o)  =  27r,  giving 


m 


The  expression  derived  for 
to  n  is 


Pi 


sin^  ^0 
sin^O 


Pi 


D 


which  is  to  be  minimized  with  respect 


R’^riy— mi 


(9.3) 


In  order  to  facilitate  the  work,  the  following  expressions  will  be  defined 


mi  = 


ni^=x 

Et 


(9.4) 


Utilizing  equation  (9.4),  equation  (9.3)  becomes 

P  sin^  ^0  ^^(a+xY  . _ ^  m  5) 

^Dsin^e  x—a  x-a 


The  expression  on  the  left-hand  side  is  to  be  minimized  with  respect 
to  X.  Therefore, 

/„  sin^0o.^l\  _(x-~a){2){x^a)~[{x-\-ay 

X 

a;^ — 2a.®  —  (to'" + a)  _ 

{x-aY 
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giving  _ 

x^a  (^1  ±  'y/ 

But  hence  the  positive  sign  must  be  chosen,  and 


x=a(l+^4+^) 


(9.6) 


To  assure  that  this  does  yield  a  minimum,  further  differentiation 
yields 


sin^ 
sin^  6 


2 


>0 


Thus,  equation  (9.6)  is  the  required  value  of  x  (or  to  obtain 
minimal  load  conditions  for  P^. 

Equation  (9.6)  may  also  be  used  to  obtain  the  number  of  buckles. 
Define  the  quantities 

'-f 

then 

7 


a__  Et  m  EtE^  ^ 

DR^  y^  EP  ~  y^ 

giving  the  buckle  number,  n,  as 

/i^=7sm20(l+2^1  +  62(^))  (9.7) 

In  order  to  substitute  equation  (9.6)  into  (9.5)  and  determine  Pi 
(min),  we  first  determine 


xpa~2a 


(‘+V'+5')"S(‘+V'+ 


6V\ 

tV 


Now  define 
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giving 

x-”a=^Vl+2/j  (l+Vl+2/)  (9-S) 


Substituting  equation  (9.8)  into  (9.5),  and  observing  that  the  maxi¬ 
mum  buckling  stresses  occur  at  6=60,  then 


p,  %  (i+Vi^)' 

Svr+-^ 


Et  12(1 
Et^ 


1 

|lVl+2/ 


But, 


2y  (1+a/1+1/)^  I  1 

VI+^  T<^Vi+y 


Pi 

D 


Pi 

Et 


12(1-.^)  = 


4^1 


Et^ 


giving 

P,_  7  P  (l+Vl+y)^  6"  P  1 

¥  VI+l/  VT+y 

7  (H~Vi+y)^  I _ L  ^ 

26  V*  Vi+y  27^14-2/ 

_j_  (i+Vi+y)T  1  1 

272/  Vl+y  27  Vl+2/ 

since 

j^_l 

y 


Now  define  ijl-py=Z,  y=Z‘—\,  giving 

Pi^JL  fi^+Zr  ,  .1  .  1  2g+2Z 

P/“27zL2^-1  J  27Z  Z^-l 


-  (9.9) 

Et  y  Z—1  7  Vl+2/— 1 

since  Z-j-l  ^0. 

Equation  (9.9)  is  brought  into  its  clearest  form  when  expanded  by 
the  binomial  theorem.  The  expansion  will  be  made  in  powers  of 
1/y  since  y  is  large. 


THE  TENSILE  PARADOX 


83 


Fi  1  1 

Et  7  1/2  /  1 

y 

1  1 
“7  r,  I V  ^  I  (1/2)  (-1/2) 

7  L  ^2  2! 


where,  for  11=0.3,  b^=3(l  —  p‘)=2.73,  and  6=1.652.  Thus, 


Pi_  1 


Et 

1.652p 

.^+^6(9 

1'  ^  { 

-Y+ . . . 

p/  j 

PiP 

'  59.5  V, 

1 

7 

1.652 

[^+5.46  (p) 

'  59.5  Vi 

■  ■  •] 

_7 

P 

giving,  finally, 


PiR 

EP 


u_ 

0.606 

1  / 

7\,  1  /7V 

1  / 

7V  , 

(9.10) 


for  the  buckling  edge  load,  Pi,  and 


71^=7  sin^ 


1+2^1+2.73  0'  (9.11) 


for  the  number  of  buckles,  n.  Since  the  buckles  should  initiate  at  the 
small  end,  6^  is  used. 

It  will  be  observed  that  7/p  increases  as  Oq  and  t  increase,  which 
means  that  larger  buckling  values  are  obtained  for  close-truncated 
thick  specimens.  In  fact,  equation  (9.9)  shows  that  Pi  is  monotoni- 
cally  increasmg,  and  approaches  infinity  for  large  values  of  7/p  since 
7/^0  when  7/p^oo  .  Such  geometries  are  approaching  that  of  a  tensile 
test  specimen  in  the  same  manner  that  short  thick  Euler  columns  are 
approaching  that  of  a  compressive  test  specimen. 

Equation  (9.10)  is  to  be  interpreted  as  a  “classical’ '  result  obtained 
from  linear  theory  and  as  not  including  the  effects  of  imperfections. 
Similar  classical  values,  using  equations  (9.1)  and  (9.2),  have  been 
obtained  for  the  sphere  under  external  pressure,  the  end  loaded  cylin¬ 
der,  and  the  cylinder  under  hydrostatic  pressure.  These  three  cases 
provide  additional  verification  of  equations  (9.1)  and  (9.2). 

Just  as  in  almost  every  other  shell  buckling  solution,  actual  design 
allowables  must  be  obtained  and  justified  by  test.  The  effects  of 
nonlinearity,  imperfections  and  test  methods,  all  of  which  can  have 
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FIGURE  9,2 — Test  values  of  C  vs  yip. 

significant  influence  upon  the  actual  test  values,  are  not  readily  in- 
corporative  in  the  theory.  This  is  the  situation  in  general  and  this 
report  is  no  exception  to  the  general  rule. 

Test  values  are  reproduced  as  a  graph  of  C  versus  yjp  in  figure  9.2. 
It  is  to  be  noted  that  more  scatter  is  to  be  expected  in  the  low  7/p 
values,  since  imperfections  have  more  influence  upon  thin  than  upon 
thick  specimens  for  a  given  radius.  Design  allowables  based  upon  a 
statistical  probability  exception  from  these  seven  tests  are  not  recom¬ 
mended  since  the  wide  scatter  in  these  few  tests  would  result  in  exces- 


THE  TENSILE  PARADOX 


85 


sively  low  C  values  and  impose  undue  penalties  upon  the  designer. 
It  is  recommended  that  further  testing  be  undertaken.  However, 
until  additional  values  can  be  obtained,  a  C  not  greater  than  the  mini¬ 
mum,  6^=0.273,  should  be  employed. 

A  final  remark  regarding  the  test  values  is  in  order.  The  buckling 
test  loads  noted  are  those  given  at  the  instant  of  the  first  buckle 
formation.  Additional  loading  results  in  mor^  complete  buckle  pat¬ 
terns  at  considerably  higher  loads. 

The  number  of  buckles  will  now  be  computed  to  verify  the  assump¬ 
tion  that  a  large  number  of  buckles  are  formed.  First,  however,  a 
close  approximation  to  equation  (9.11)  will  be  derived.  Since 

2.73  (^y»l 

we  may  write 

7^2=7  sin^  ^0  j^l+2V^3 

If,  in  addition, 

2V^  (^)»1 

we  obtain 

7i=1.82Vp  sin  ^0  (9.12) 


TABLE  9.1 — Approximate  n  Values  From  Calculation  and  Test 


Spec. 

7 

p 

p 

7 

From 

(9.10) 

Value  of 

From 

(9.9) 

n 

Approx., 
as  meas- 
sured  » 

1 

57.6 

481 

8.35 

36.6 

37.3 

30 

2 

57.6 

477 

8.  27 

36.4 

37.  2 

30 

3 

57.6 

454 

7.88 

35.  6 

35.7 

30 

4 

27.3 

759 

27.8 

41.3 

41.6 

28 

5 

27.3 

638 

23.4 

37.9 

38.2 

30 

6 

57.6 

1600 

27.8 

66.8 

67.2 

50 

7 

57.6 

1510 

26.  2 

64.8 

65.3 

35 

a  Only  approximate  n  values  as  measured  from  tests  could  be  obtained, 
since  complete  buckle  patterns  were  not  formed  in  some  cases,  and  in 
others,  estimates  had  to  be  made  from  photographs. 
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